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1 Introduction 



C^l ' The whole set of correlation functions fully describes a Quantum Field Theory, as it is related to the S- 

matrix elements. In QCD or pure Yang-Mills theories Green functions are most often gauge dependent 
quantities which have no direct relationship with physical observables, the latter being necessarily 
gauge invariant. However, their indirect physical relevance is well known. In particular, long distance 
(or small momentum) Green functions will hopefully shed some light on the deepest mysteries of QCD 
such as confinement, spontaneous chiral symmetry breaking, etc. 

Indeed, more than thirty five years after the discovery of QCD and notwithstanding its numerous 
successes, notwithstanding cither the fact that everybody is convinced that QCD implies confinement, 
a real proof of it from first principles has not yet been achieved. This is doubtless one of the major 
scientific challenges of this century. In fact, in the case of QCD, one is not even provided with a 
precise mathematical formulation of confinement. In the case of the pure Yang-Mills theory, such a 
mathematical formulation at least exists: it is the area-law for Wilson loops which unluckily does not 
hold for QCD due to the string breaking by the sea quarks. Furthermore, the Polyakov loop (the product 
of link variables along a curve in time direction closed by periodic boundary conditions) , probing the 
screening properties of a static colour triplet test charge, appears to be the order parameter for the 
deconfinement phase transition. Here we will restrict ourselves to the pure Yang-Mills theory. We will 
not enumerate all the tracks which have been followed to understand confinement by means of the 
peculiarities of QCD in the infrared domain, there are many of them (see, for instance, the classical 
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Wilson's |l| or Cornwall's |2j works, or the very recent introduction to the confinement problem in 
ref. Q and references therein). In this review we will only mention the Kugo-Ojima and the Gribov- 
Zwanziger approaches. The latter aimed in principle to deal with the problem of Gribov copies, but 
it is also thought to be connected with confinement scenarios. Both approaches arc related to the 
behaviour of the Green functions of the gluons and ghosts in the deep infrared. 

It is now well established that the vacuum of a quantum field theory is never trivial, and especially 
not in the case of a non-Abclian gauge theory. It is believed that the low modes of the vacuum, 
the condensates, are the keys to understand its non-perturbative properties. The vacuum is a gauge 
invariant state. But the configurations of the fields in the vacuum have different features in different 
gauges. This allows for different, complementary and rewarding views into its properties. We will 
address this issue (sec. I3.8|) . 

We will concentrate our efforts on a review of the gluon and ghost Green functions properties 
at small momentum in a pure Yang-Mills theory We will invoke results from the analysis of Dyson- 
Schwinger equations (DSE), Slavnov- Taylor identities (ST) and from lattice QCD, paying special at- 
tention to the interplay of all these techniques and particularly to the study of Yang-Mills solutions 
when lattice results happen to be applied as DSE inputs. 

We work hereafter mainly in the Landau gauge, but some of the results we will present are actually 
valid in any covariant gauge. We will also discuss some results in Coulomb gauge and focus especially 
on the similarities of the general properties of the solutions in both Landau and Coulomb gauges. 



Our notations will be the following : 

(F^) ab (k) = S ab ^p- (1) 

{G ^ (k) = 5ab ^^ v _b^j (2) 

r$ c p (p,q,r) = f abc r^ p (p,q,r) (3) 

r c ; bc (q 7 k;p) = r bc (- lqv )g r^(q,k;p) 

= mf abc ( -q v H x {-q,k)+p v H 2 {-q,k) ) , (4) 



respectively for the ghost propagator, the gluon propagator, the three-gluons vertex and the ghost- 
gluon vertcjIH. All momenta are taken as entering. In cq. (U) —q is the momentum of the outgoing ghost, 
k the momentum of the incoming one and p = —q — k the momentum of the gluon. i?i(2) corresponds 
to the gluon-transverse (longitudinal) form factor that will be extensively invoked in the following. g$ 
is the bare coupling that should be properly rcnormalized and, in the limit of vanishing incoming 
ghost momentum (k — > 0), one has 

M<7 2 ) = ^F 2 (q 2 )G(q 2 ) , (5) 

which gives the running coupling in the Taylor renormalization scheme, where F(p 2 ) and G(p 2 ), de- 
fined by Eas. (|ll2[) . are the dressing functions of the ghost and gluon propagators respectively. Up to 
logarithms, we parametrise the propagators in the infrared by setting at leading order 

CtG 

(6) 

F(p 2 ) — I — ) , when p 2 is small, 

where A, rj are some dimensional parameters. Finally we shall set D(q 2 ) = G{q 2 )/q 2 

As we mentioned, our goal will be to describe the current state-of-the-art concerning the low- 
momentum properties of the Green functions in gluo-dynamics, mainly by focussing on the results 

1 We stick to the decomposition given in ref. [4| except for the arguments of the scalar functions, for which 
we keep the same order as in r itself 

2 By Z'g 1 — Zi 1 , where Zz{Zz) is the gluon (ghost) propagator renormalization constant and Z\ is 
the proper ghost-gluon vertex renormalization constant. 
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obtained with two main approaches: DSEs and Lattice QCD. The paper is organized as follows: a first 
section (sec. [5]) will be devoted to provide the reader with a first insight on the problem of the Yang- 
Mills Green functions low-momentum behaviour through reviewing the main analytical properties 
recently derived in the literature for the solutions; we will review then, as exhaustively as possible, 
the numerical lattice QCD results on the subject in sec. [3j the results obtained from the numerical 
resolution of DSEs, with different truncation approaches, will be discussed and put in connection with 
lattice results in sec. @] and we finally conclude in sec. [5j However, we will start now by "saying" a few 
words to introduce those two main approaches and their results, but also about other approaches as 
the ones based on redefining the QCD lagrangian to properly correct the Gribov ambiguity. 

1.1 Lattice results 

The lattice technique in gauge field theories was initiated by Wilson in 1974. It has now been used for 
the field of Yang Mills Green's functions for nearly 25 years. Let us begin by recalling what it consists 
in. The general idea is to rewrite the lagrangian on a discretiscd space-time (the lattice) instead of the 
ordinary continuous one. Accordingly all derivatives (including the covariant derivatives) are replaced 
by finite differences. If in addition one limits oneself to a finite volume, the path integral is now over 
a finite (but huge) number of field variables and can be estimated as a sum over a stochastically 
determined set of configurations. The advantages are well known : 

• The only ingredient is the original lagrangian, excluding any additional hypothesis. So, in some 
sense, the results obtained in this way can be considered as exact from the theoretical point of 
view. 

• The technique is essentially non-perturbative. 

• The lattice spacing a acts as an ultraviolet cutoff : no divergence occurs. The continuum limit can 
be recovered by letting a go to zero, at the price of an appropriate renormalisation. 

There are some drawbacks to the method, however. 

• It is very demanding in computing power. 

• The results are given in numerical form and, as such, they suffer from uncertainties. The first one 
is of statistical nature and stems from the finite number of configurations which are used in the 
evaluation ; it can in principle be reduced at will by increasing this number. The other sources are 
more intrisically related to the technique, they give rise to systematic errors. We mention them 
briefly below although they will be discussed in greater detail in dedicated sections in the following 
(see section ET5"|) . 

First of all the actual calculations are necessarily performed at finite values of the lattice spacing. 
Naively, at tree level, the discretisation procedure generates effects of order a in the fcrmionic lagrangian 
and of order a 2 in the Yang-Mills one. The gap between those conditions and the continuum limit can 
be partly reduced by using an improved lagrangian (see for instance refs. [|[ for the fermionic case and 
@ @| i n the pure gauge one) . The remaining part has to be treated numerically. 

Second, going to the lattice implies, as we have already said, that one works in a finite volume. This 
induces, in turn, a 1/L spacing in p-space. Since we are interested in the infrared limit of the Green's 
functions we shall have to take special care of potential discontinuities in the neighbourhood of 0. 

Third, the space on which we are working is actually a torus since the field configurations are usually 
chosen to be periodic. One can therefore wonder what kind of relationship the functions obtained in 
this way entertain with the "real world" ones. 

As a final remark, we note that the Gribov copies problem, which is absolutely general, can be 
specifically dealt with on the lattice. The gauge fixing procedure consists in gauge-transforming the 
field configuration into one of its gauge equivalents which satisfy the given gauge condition. Generally 
the different possible choices do not contribute equally to the path integral. It is therefore an important 
issue to know how large this effect is and how one can minimize it. Within the numerical precision, 
this can be explicitely done on the lattice, as will be discussed in section 13.61 

The number of papers dealing with lattice simulations for Green functions in the IR is rather 
large, and their results are given according to quite a variety of presentations, which makes the task 
of comparing them not so easy. We shall try to give an account of those works in sec. |3] Note that the 
most recent studies have been performed in the unquenched theory (0; Qjl El; El), so that they fall 
out of the scope of this review and will not be discussed. 
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1.2 The DSE picture 

In principle, the field equations for the Green functions of a quantum field theory can be derived from 
the integral representation of the theory; these are the Dyson-Schwinger equations (DSEs) of this theory 
which are considered to describe its non-perturbative dynamics fL3t [l4| . Under the only assumption 
of the existence of a well-defined measure in a functional integral representation of the generating 
functional for the Green functions of a theory, QCD for instance, the corresponding DSEs can be 
derived (see for instance ^3) as an infinite tower of coupled integral equations. Of course, in practice, 
these equations need to be truncated to be studied. Typically, one can make appeal to additional 
sources of information, for instance Slavnov- Taylor identities, to express higher n-point functions in 
terms of elementary two-point ones or to generate some general ansatze for vertices. Thus, a closed 
coupled system of DSEs is to be obtained and can be numerically solved, subject to the validity of the 
truncation rules which are applied. This DSE approach has been extensively applied to investigate the 
low-momentum behaviour of the QCD Green functions; in particular, for the Yang-Mills gluon and 
ghost propagators. We will dedicate sees. [5] and @] to describe some of the main results contributing 
to the current DSE picture for this low-momentum behaviour of Yang-Mills Green functions. Let us 
however briefly introduce those results and how they emerged in the last few years. 

After the first pioneering works of Mandelstam [lfl [l7| and the further ones of Brown and Penning- 



ton 18(, only a few years ago, a new paradigm emerged when it was widely accepted (see for instance 
[TEl ]) that a vanishing gluon propagator and a diverging ghost dressing function at zero- momentum in 
Landau gauge made up the unique solution of the truncated tower of DSEs. In contrast, alternative 
DSE solutions were also predicted to give a massive gluon propagator [20|. Lattice QCD (LQCD) 
estimates for those propagators appeared to be also in contradiction with a gluon propagator that 
vanishes at zero- momentum or with a ghost dressing function that diverges [2ll : l22l l23l: [24 : [25j. We 
addressed this issue in two recent papers [2(1 : l27l | and tried to clarify the contradiction. After assuming 
in the vanishing momentum limit a ghost dressing function behaving as F(q 2 ) ~ [q 2 ) ap and a gluon 
propagator as D(q 2 ) ~ (g 2 ) QG_1 (or, by following a notation commonly used, a gluon dressing function 
as G(q 2 ) = q 2 D(q 2 ) ~ (q 2 ) aa ), we proved that the ghost propagator DSE (GPDSE) admits two types 
of solutions: 

• If cxf 7^ 0, the low-momentum behaviour of both gluon and ghost propagators are related by the 
condition 2o.f + Qg = implying that F 2 {q 2 )G(q 2 ) goes to a non- vanishing constant when q 2 0. 
This solution is now called "scaling" . 

• If ctF = 0, the low-momentum leading term of the gluon propagator is not constrained any longer 
by the leading but instead by the next-to-leading one of the ghost propagator, and LQCD solutions 
indicating that F 2 (q 2 )G(q 2 ) — > when q 2 — > [23|; [24[ can be pretty well accommodated within 
this case. This solutions is named as "decoupling . 

In particular, the numerical study in ref. (2(| of the GPDSE using a LQCD gluon input finds that 
two classes of solutions emerge, depending on the value of the strong coupling constant at the renormal- 
ization point, which is a free parameter in this exercise. Indeed, it seems to be by now well established 
that the two classes of solutions, decoupling and scaling may emerge from the tower of DSE [l9l: I20I; l28l| . 
Such a nomenclature, despite being widely accepted, can be misleading. The perturbativc running for 
the coupling constant renormalized in Taylor-scheme is given by Eq. ([5]) and one can thus extend this 
definition, although not univocally, to the IR domain. However, a scale invariance ( a decoupling) of the 
IR dynamics for the theory cannot be inferred from the low-momentum behaviour of such a coupling 
in the scaling (decoupling) solution. In particular, as will be seen in next Eq. (joT)]) of subsection 12.21 
(see also Eqs. (2.30-2.32) in ref. [29j). an effective charge can be properly defined for phcnomcnological 
purposes such that it reaches a constant at zero-momentum in the "decoupling" case, which means the 
absence of decoupling. Leaving aside nomenclature, the two classes of solutions are definitely different 
and the Taylor-scheme coupling, although maybe not appropriate for phcnomcnological purposes in 
the IR domain, is a well defined and very convenient quantity to discriminate them. 

How both types of IR solutions for Landau gauge DSE emerge and how the transition between 
them occurs, in relation with the size of the coupling (taken as an integration boundary condition at 
the renormalization momentum), was initially discussed in ref. [26[ through the analysis of a ghost 
propagator DSE combined with a gluon propagator taken from lattice computations. It should be 
remembered that one needs to know the QCD mass scale to predict the QCD coupling at any mo- 
mentum. This mass scale should be of course supplied to get a particular solution from DSE and can 
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be univocally related to the boundary condition needed, after applying a truncation scheme, to solve 
the equation. The existence of a critical value for the coupling at any renormalization momentum was 
suggested by that partial analysis. No solution was proven to exist for any coupling bigger than the 
critical one and the unique scaling solution^] seemed to emerge when the coupling took this critical 
value. Later, the authors of ref. [28[ confirmed, by the analysis of the tower of DSE truncated within 
two different schemes and also in the framework of the functional renormalization group, that the 
boundary condition for the DSE integration determined whether a decoupling or the scaling solution 
occurs. A similar analysis has been recently done in the Coulomb gauge [3l| leading to the same pat- 
tern as in ref. [2(| , although the authors interpret the boundary condition in terms of the gauge- fixing 
ambiguity (see also [H|;[33j]). Furthermore, an analytic study based on the pinch technique (PT) in 
ref. [341 ] shows that, within some approximations, there is a lower limit for the gluon mass, below which 
the PT coupling is singular in the IR, that can be also interpreted as an upper limit to the coupling at 
some renormalization point. Very recently also, a next-to-leading low-momentum asymptotic formula 
for the decoupling ghost dressing function solutions was obtained by studying the ghost propagator 
DSE under the assumption, for the truncation, of a constant ghost-gluon vertex and of a simple model 
for a massive gluon propagator [35| . In this asymptotic formula, the ghost-propagator low- momentum 
behaviour appears to be regulated by the zero-momentum effective charge in Taylor scheme [29l | and 
by the Landau-gauge gluon mass scale. 

That DSEs, being an intricated tower of coupled functional differential equations, admit multiple 
solutions belonging to different types is not very surprising. The curent lattice results appear to be 
clearly compatible with only one type of solutions and, after the appropriate physical calibration of 
the simulations has been performed (and all the lattice artefact have been put properly under control), 
must help to select the "-physical 1 solutions for the QCD Green functions among the multiple DSE 
ones. Furthermore, as will be discussed in the next subsection, an effective action incorporating the 
so-called Gribov horizon can be properly built to deal with the Gribov problem, and the lattice-like 
solutions can also appear as those minimizing this effective action (the emergence of non-vanishing 
dimension-two condensates, that we will also discuss in section 13.81 plays a crucial role in obtaining 
the lattice-like solutions). 

1.3 The Gribov and the Gribov-Zwanziger approaches 

The computation of gauge dependent quantities like the correlation functions we are interested in 
suffers from the Gribov ambiguity: imposing a constraint like the Landau gauge condition d^A 11 = 
is not enough to pick unambiguously a unique representative from each gauge orbit (36| . Gribov then 
proposed to reinforce the gauge condition by restricting the space of gauge fields to representatives 
minimizing the functional 

T(<P)=Tr J dxA*(x)A*(x) (7) 

where A* is the gauge transformed field. This is appealing from two points of view : first the Gri- 
bov problem is reduced as gauge copies which are maxima or saddle points of T are pushed out of 
the game and second this restricted space, named the Gribov region, is rather easy to characterize 
mathematically. Actually the second derivatives of the functional T defines an operator, the so-called 
Fadcev-Popov operator, which has to be positive in the Gribov region. For the Landau gauge this 
operator reads : 

M ab = QHjjab ) (g) 

where D is the covariant derivative. But in general along a gauge orbit the minimum for the functional 
is not unique. The ambiguity in the gauge fixing is not fully eliminated this way. It is then natural to 
try to restrict further the space of gauge fields to configurations which are not only a minimum but an 
absolute minimum of J- . This region, named the fundamental modular region, in which on each gauge 
orbit F possesses a unique absolute minimum (cf reference [13]), would be theoretically convenient to 

3 The authors of [3(| proved there, once the scaling behaviour is assumed, the uniqueness for Yang-Mills 
infrared solutions 
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eliminate the Gribov problem but unfortunately no practical explicit characterization has been found 
so far. 

Gribov proposed a heuristic method to perform the restriction to the Gribov region. He computed 
the inverse of the Faddeev-Popov operator in perturbation theory at one loop and found an expression 
with a structure like : 

1 ' fc 2 l-a(k,A) [ ' 

The explicit form for a can be found in the original paper by Gribov but is not necessary for the general 
discussion here. The Gribov region is a space where the eigenvalues of the Fadeev-Popov operator are 
positive. Forbidding the crossing of the boundary of this region is equivalent to forbid the appearance 
of a zero eigenvalue. This is the origin of a condition proposed by Gribov to restrict the space of gauge 
fields to the Gribov region : 

a(0,A) < 1 (10) 

With this restriction as a constraint, Gribov computed the ghost and the gluon propagators (at one 
loop) and found : 

k 4 

G Grlb ov(k 2 ) = fc4 + to 4 (11) 

FGribov(k ) - (12) 

Namely, when k 2 goes to 0, a vanishing gluon propagator and an enhanced ghost propagator are 
exhibited, vtiq is called the Gribov mass. 

One step further has been accomplished by Zwanziger [38| to relax the approximation used in the 
determination of the constraint and to extend it to all orders. He found a condition to restrict the 
space of gauge fields [38| 

h(A) < d(N^-l) (13) 

where h(A), the horizon function , is a non local functional of the gauge fields. This approach has been 
subsequently refined and developed by Zwanziger himself ([39||) and other authors ([1(1 5J; H3; IH|). 
Zwanziger showed that this condition can be exponentiated to be transformed into a contribution Sgz 
to the action : 

S GZ = 7 (h(A)-d(N c -l)) (14) 

to be added to the usual contributions coming from the gauge action and the gauge fixing terms. 7 is 
defined by a relation, the horizon condition, required to implement the exponentiation : 

< h(A) > GZ = d(N 2 - 1) (15) 

This action has many interesting properties : auxiliary fields can be introduced to transform Sgz in a 
local form [44[ ; the theory has been proven to be renormalizable [44| ; the gluon and ghost propagators 
have been computed and results similar to the Gribov ones are found: a vanishing gluon propagator and 
an enhanced ghost propagator. These results were in accordance with the common prejudice inferred 
at that time from the Dyson-Schwinger equations. 

As it became clearer and clearer that the lattice results for the propagators were not in accordance 
with these predicted behaviors, a refined version of the GZ approach was necessary to have an answer 
for this discrepancy. The key point is the explicit introduction of the dimension two condensates [4ll:l43j. 
This has the very nice feature that the renormalizability is not spoiled. With this refinement the gluon 
and ghost propagators are in qualitative agreement with the lattice results and the prediction from 
the "decoupling" solution of the Dyson-Schwinger equations. In particular, the gluon propagator in 
the so-called "refined" Gribov-Zwanziger (RGZ) formalism is shown to behave as 

2 k 2 + M 2 

D(yk ' ~ k 4 + k 2 (m 2 + M 2 ) + 2g 2 N c l 2 + M 2 m 2 ^ 
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where 7 is the Gribov-Zwanziger parameter determined by the horizon condition, Eq. (|15|) . while M 
and m are two mass parameters related to the two above mentioned condensates . In particular, m 
is related to the dimension- two gluon condensate, (A 2 ), that will be discussed in detail below (cf. 
section I3.8|) and plays a crucial role for the RGZ gluon propagator to account for lattice results [45( . 
Furthermore, the authors of ref. [46^ studied the effective action in RGZ and provided firm evidences 
that the dimension-two condensates should be non-vanishing (and hence the mass parameters M and 
m) and a strong indication that, in addition to the non-zero gluon propagator at vanishing momentum, 
the ghost propagator is not enhanced, in consistence with lattice data and with the previously discussed 
DSE picture for decoupling solutions. 



1.4 Other approaches 

Many other approaches have been applied to investigate the low-momentum properties of the Yang- 
Mills Green functions. We will end this introduction by indicating some of them and addressing the 
interested reader to some original works. Apart from DSEs or Lattice QCD, one of the most followed 
approaches is that of functional renormalization group equations (FRGs) [47| HH |49|; HO; [5l|. Indeed, 
FRGs and DSEs appear to be rather interconnected and, for instance, it has been shown that the 
integrated flow equations define a set of DSEs within a particular renormalization scheme [5lj . First, 
it was reported that FRGs analysis of the Landau-gauge low-momentum behaviour for the Yang-Mills 
Green functions leads to the uniqueness of a scaling-type solution [5(1 H(| , but it has been recently 
proven, also within FRGs, that both scaling and decoupling solutions exist (52[. On the other hand, 
other approaches like the infrared mapping of A0 4 and Yang-Mills theories in ref. [5j| H3; HE1 or 
the massive extension of the Fadeev-Popov action in ref. [56l ; \El\ appear to support a massive gluon 
propagator and a free ghost, i.e. a decoupling-type solution. In particular, a very accurate description 
of lattice data for Yang-Mills gluon and ghost correlators in Landau gauge is obtained by means of a 
one-loop computation with this last massive extension of the Fadeev-Popov action in four dimensions, 
while the main features for lattice results in d=2,3 can be also accounted for. Furthermore, the low- 
momentum behaviour for the one-loop results with this massive action matched pretty well with the 
one expected for a decoupling solution within DSEs approach. One can find in the literature still other 
methods, such as the application of stochastic quantization [58|; HH, that have been applied to the 
subject. 



1.5 The low-momentum correlators and the confinement problem 

Let us end this introduction with one comment about the Kugo-Ojima (60| confinement criterion which 
might be in order here. In the Kugo-Ojima colour confinement picture, the physical spectrum of the 
theory is free of coloured asymptotic states as a consequence of the so-called "quartet mechanism" . 
A sufficient condition for it to take place is that a certain correlation function, usually denoted by 
u(q 2 ) and called "Kugo function", should satisfy u(0) = — 1. Furthermore, in Landau gauge, the Kugo 
function is linked to the ghost dressing function such that F(0)(1 + u(0)) = 1 (This relation was first 
noted by T. Kugo [6l| and very recentely, K-I. Kondo triggered an interesting discussion about this 
relation, in connection with the Gribov horizon condition and its implications on the Landau-gauge 
Yang-Mills infrared solutions |4C| |62|). Therefore, the sufficient condition for the realization of Kugo- 
Ojima confinement scenario requires a divergent ghost dressing function. Then, for the two classes of 
solutions above discussed, only the scaling may satisfy this sufficient condition. A different mechanism 
should thus explain the confinement for decoupling solutions, as for instance the one provided by the 
center vortices scenario @; |&| H3] • 

However, as stated in the brief review about "Strong Coupling Continuum QCD" recentely appeared 
in the proceedings of the 9th conference on "Quark Confinement and the Hadron Spectrum" [65|], the 
physics of hadrons and the confinement do not depend very much on how gluon and ghost propagates 
over distances of atomic scales but mainly on those of the size of a nucleus, where the difference between 
scaling and decoupling solutions is not important. Thus, it is very reasonable to think that the real 
QCD confinement mechanism might not be of a great help to discriminate among the two types of 
low-momentum solutions. 
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2 A first insight with analytical tools 

2.1 The low momentum solutions from the ghost propagator DSE 
2.1.1 The ghost propagator DSE 

We will examine the Dyson-Schwinger equation for the ghost propagator (GPDSE) which can be 
written diagrammatically as 

q-k 



o 

k 



V 




a,k c,q 



i.e., denoting by (resp. G^) the full ghost (resp. gluon) propagator, 



(F (2) )a b 1 (k) = -S ab k 2 

9ofacdfebj 



d 4 q 



(q) (iq v , )r u , u (-q,k;q-k)^)(G^)li(q-k), 



where F stands for the bare ghost-gluon vertex (above defined in Eq. ((][])), 

rf c (-q 7 k;q-k) = zg a f abc q v ,f u , v {-q,k ]q -k) 

= ig a f abc (q l/ H 1 (q 1 k) + (q-k) l/ H 2 (q 1 k) ) , 



(17) 



(18) 



while q and k are respectively the outgoing and incoming ghost momenta and go is the bare coupling 
constant. Let us now consider eq. (fTT)) at small momenta k. After applying the decomposition for the 
ghost-gluon vertex in eq. (|18[) . omitting colour indices and dividing both sides by fc 2 , it reads 



1 



F(fc 2 



, d A q [ F(q*)G((q-k?) 

i + g N c 



(k ■ qf 



Hi(q,k) 



(19) 



It should be noticed that, because of the transversality condition, H2 defined in eq. (JT5J) docs not 
contribute for the GPDSE in the Landau gauge. 



2.1.2 Renormalization of the Dyson-Schwinger equation 

The integral equation eq. (|19p is written in terms of bare Green functions. It is actually meaningless 
unless one specifies some appropriate UV-cutoff0 A, and performs the replacements F{k 2 ) — > F(k 2 , A) 
. . . .It can be cast into a renormalized form by dealing properly with UV divergences, i.e. 

9U^)^Z- 2 ( t , 2 ,A)g 2 (A) 



G R (k 2 , f i 2 )^Z~ 1 ( fJ ; 2 ,A)G(k 2 ,A) 
F R (k 2 }f i 2 )=,Z^^ 2 7 A)F(k 2 ,A) , 



(20) 



where p? is the renormalization momentum and Z g) Z^ and Z% the renormalization constants for 
the coupling constant, the gluon and the ghost respectively. Z g is related to the ghost-gluon vertex 

renormalization constant (defined by Fr = Z\Fb) through Z g = Zi(Z^ 2 Z^)~ x . Then Taylor's non- 
renormalization theorem, which states that Hi(q, 0) + H2(q, 0) = 1 in Landau gauge (see app. lA"| and to 

4 We have written for simplicity the UV cutoff as a hard cut-off. It is preferable to use a gauge inv arian t 
regularization procedure in view of the advantage of exploiting Ward-Slavnov- Taylor identities (see sec. 12.31) . 
In practice we will derive our results from the subtracted GPDSE which incorporates gauge invariant UV 
regularisation. 
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any perturbative order, can be invoked to conclude that Z\ is finite. We recall that the renormalization 
point is arbitrary, except for the special value u = which cannot be chosen without a loss of generality 
(see, in this respect, the discussion in ref. j66|)- Thus, 



1 



F R {k 2 ^ 2 ) 



Z 3 (li 2 ,A)+N c Z 1 g 2 R (^ 2 ) E R {k 2 ^ 2 -A) 



(21) 



where 



S R (k 2 ,fi 2 ;A) 



d 4 q 



(27T) 4 

F R {q 2 ^ 2 )G R {{q-k) 2 ^ 2 ) 



(fc ■ qf 
k 2 



H 1:R (q,k;n 2 ) 



(22) 



One should notice that the UV cut-off, A, is still required as an upper integration bound in cq. (1221) 
since the integral is UV-divergent, behaving as J dq 2 /q 2 (l + llas/(27r) log {q/ ^)))~ 35 ^ 44 ■ In fact, the 
cut-off dependence this induces in S R cancels against the one of Z 3 in the r.h.s. of eq. (|2"Tj) . in 
accordance with the fact that the l.h.s. does not depend on A. 

Now, we will apply a MOM renormalization prescription. This means that all the Green functions 
take their tree-level value at the renormalization point and thus: 



F R (n 2 ,V 2 



G R (^ 2 ^ 2 



(23) 



In the following, H\{q,k) will be approximated by a constant with respect to both momenta and, 
provided that Hi(q, 0) = I at tree-level, our MOM prescription implies that H± jR (k, q; /i 2 ) = 1 and Z\ 
is a constant in terms of fi. 



2.1.3 A subtracted Dyson-Schwinger equation 



The renormalized GPDSE, eq. (f2~Tj) . should be carefully analysed. We aim to study the infrared be- 
haviour of its solutions and therefore focus our analysis on the momentum region, k <C ^IqcD; where 
the IR behaviour of the dressing functions (presumably in powers of the momentum) is supposed to 
hold. One cannot forget, though, that the UV cut-off dependences in both sides of eq. (f2~T|) match 
only in virtue of the previously mentioned relation between the ghost and gluon propagator anomalous 
dimension and the beta function. 

However, in order not to have to deal with the UV cut-off, we prefer to approach the study of the 
GPDSE in the following manner: we consider eq. (|2ip for two different scales, Afc and \nk (with k < I 
some fixed number and A an extra parameter that we shall ultimately let go to 0) and subtract them 



F R (X 2 k 2 ,fi 2 ) F R (X 2 n 2 k 2 ,n 2 ) 



N c g R {f) Zi I A7 fl (A 2 fc 2 ,^;co) - S R (X Z K z k' , ^ ; ^) 



2 „2 .2 



(24) 



5 One can easily check that Z 3 1 (/i 2 , A)En(k 2 , p 2 ; A) approaches some finite limit as A — > oo since the ghost 
and gluon propagator anomalous dimensions and the beta function verify the relation 27 + 7 + /3 = |67l |. 

6 This approximation is very usually used to solve GPDSE. Some lattice data are available for the gh ost-g luon 
vertex, although they are by far less numerous than the ones regarding the propagators (see section The 
pr ese nt data do indicate that the zero gluon momentum H\{q,q) is approximatively constant with respect to 
q |23l j. Of course, more data for different kinematical configurations should be welcome to check that approxi- 
mation. 
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Then the integral in the r.h.s. is UV-safe, thanks to the subtraction, and the limit A — > oo can be 
cxplicitely taken, 



£ R {rk 2 , oo) - £r(\ k k , //; oo) 



A / F{q\l?) f(k-q) 2 2 
13 1 



(27T) 



G((g-Afc) 2 , M 2 ) 
{q-Xkf 



(A Ak) 



(25) 



This equation is evidently a necessary consequence of the original one (|2Tj) . That, conversely, it is 
actually sufficient was shown in pt| . For an accurate analysis of cq. (|24|) it is convenient, in addition, 
to split the integration domain of eq. (|2"5j) into two pieces by introducing some new scale q 2 (qo, 
typically of the order of Aqcd, is a momentum scale below which the deep IR power behaviour is a 
good approximation), 



E R (X 2 k 2 ^ 2 ;^)-E R (X 2 K 2 k 2 ^ 2 ;^) = J IR (A) + / UV (A) 



(26) 



where Im represents the integral in cq. (j25|) over q 2 < q 2 and iuv over q 2 > q 2 . Only the dependence 
on A is written explicitly because we shall let it go to zero with k, n and kept fixed. The relevance 
of the q 2 scale stems from the drastic difference between the IR and UV behaviours of the integrand. 
In particular, for (Afc) 2 <C q 2 ,, the following infrared power laws, 

F m {q 2 ^ 2 ) = A^ 2 )(q 2 ) aF 

G m ((q-Xk) 2 ^ 2 ) = B^ 2 )((q-Xk) 2 ) aG , (27) 

will be applied for both dressing functions in Ttr. 

Now, a straightforward power-counting argument shows that lis, is infrared convergent if : 

ap > — 2 IR convergence at q 2 = 

a G > IR convergence at (q - fc) 2 = and (q - nk) 2 = (28) 

We shall suppose in the following that these conditions are verified. Let us first consider Iuv- Its 
dependence on A, which is explicit in the factor inside the square bracket of eq. (|2"5)) . should clearly be 
even in A : any odd power of A would imply an odd power of q ■ k whose angular integral is zero. Since 
the integrand is identically zero at A ~ and the integral is ultraviolet convergent, it is proportional to 
A 2 (unless some accidental cancellation forces it to behave as an even higher power of A) . On the other 
hand, after performing the change of variable q — > Xq, the IR contribution of the integral in eq. (|25|) 's 
r.h.s can be rewritten as: 



I m (X) ^ (x 2 ) {aF+aG) 

X \(q-k) 2 



A(v 2 )B(p 

,«G-2 



q 



2\a f 



((q - Kk) 



(27T) 

2\ QG-2 



4 («") 



(fc ■ qf 

k 2 



that, as it shall be seen in the next subsection, asymptotically behaves as 

A 2( QG +a F ) j f ac + ap < 1 

^ir(A) ~ \ X 2 In A if olq + ap = 1 



(29) 



(30) 



A J 



if qg + ap > 1 



Thus, in all the cases, the leading behaviour of Jir + iuv, as A vanishes, is given by Iir in eq. ([30 
The subtracted renormalised GPDSE reads for ac + ap < 1 as: 



1 



1 



F R {X 2 k 2 ,^ 2 ) F R (X 2 n 2 k 2 ,n 2 ) 



N c g 2 R (v 2 ) Zi Im(X) 



(31) 



for small A. We have assumed that Hi is constant when varying all the momenta but (|30I31[) remain 
true if one only assumes that Hi behaves "regularly" for q 2 , k 2 < q 2 (i.e. is free of singularities or, at 
least, of any singularity worse than logarithmic). 
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2.1-4 The integral for the ghost self-energy 

The present section is devoted to the quantitative analysis of the integral Iir(A), defined in Eq. (|29|) . 
which gives the contribution of the ghost loop to the renormalised GPDSE Eq. (|31[) . If a F + a G < 1 
it is possible to perform analytically the integral and to find a compact expression for it. In this case, 
one can write 

J IR (A) ~ A{^)B{^ 2 ) ( A 2 ) (QF+QG) ($(k; a F , a G ) - <2>( K fc; a F , a G )) (32) 
where A(/i 2 ) and B(/i 2 ) were defined in Eq. (|27p and 

#(*;«F,ao) = / -J^ ((9-fc) 2 ) QG_2 (^-9 2 ) , 03) 

provided that <P(k; ap, olg) is n °t singular, so that the subtraction inside the bracket and the integral 
operator in Eq. (|32[) commute with each other. Then, following [68j . we define 



_ r{2 + a)r(2 + b)r(-a -b-2) 

r{~a)r{-b)r{4 + a + b) ' 



(34) 



and obtain 



where 



$(k\a F ,a.G) = — ^0^2 — <t>{uF,a G ) (35) 



{a F ,a G ) = -- (f{a F ,a G - 2) + f(a F ,a G - 1) + f(a F -l,a G - 1)) 

+ \ (f(a F - 1, olg ~ 2) + /(of - 1, o G ) + /(ajr + 1, a G - 2)) . (36) 



Thus, if a F + a G < 1, 

/ir(A) * } (A 2 fc 2 )° f + " G (1 - ^») <Ko F ,o G ) . (37) 

We will now compute the leading asymptotic behavior of /ir as A — > when a F + o G = 1- In that 
case, after performing in Eq. (|29[) the following expansion, 

[(* - <z) 2 P~ 2 - [(«* - 1)T G ~ 2 - (? 2 ) aG " 2 («g - 2)(1 - «) (38) 

q ■ k , , f k 2 , . (q ■ k) 2 

-2%- + (1 + k) ^ + 2(a G -3) U 



9 2 \ q 2 9 4 

and neglecting the term odd in — > — one finds for the leading contribution 

MA) c -fc 2 (i - ^ 2 ) 2A( ( 2 y 2) A2 / 90A ^ g 2(QF+QG) - 3 

<20 sin 4 6» ^a G - 2 + 2(a G - 3)(a G - 2)cos 2 

~ fc 2 (i - K 2 ^(M 2 )y) ac(ac _ 2)A 2 ln A (39) 

3Z7T Z 

We do not specify the lower bound of the integral over q in Eq. (|3T))) because it necessarily contributes 
as a subleading term, once the ghost-loop integral is required to be IR safe. Then, as was indicated 
in anticipation in Eq. (f50|) . 7tr diverges logarithmically as A goes to zero if a F + a G = 1. In fact, 
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since Eq. (|37|) is a reliable result for any up + ac < 1 however close it may be to 1, such a divergence 
appears as a pole of a T function of 4>{otp, ag) in Eq. (1351) . 

Finally, if olf + ag > 1, the leading contribution for Iir(A) as A vanishes can be computed after 
performing back the change of integration variable, q — > g/A, in Eq. (f2U)). The first even term in Eq. (|38l) 
dominates again the expansion after integration, but now it does not diverge. Then, if we procceed as 
we did in Eq. (j39|) . we obtain 

MA) - - QG( " G = 2 \ ^"TT" 1 A{f)B{f) fc 2 A 2 (l - « 2 ) , (40) 

for small A and olq + ctp > 1. It should be noticed that /ir in Eq. (|40p depends on the additional scale 
<jo introduced in Eq. (|26[) to separate IR and UV integration domains. In fact, if one takes go — ^ OO; ^ir 
diverges. This means that, when a?i? + «g > U the behaviour of the IR power laws hampers their use 
for all momenta in the integral. The finiteness of the ghost-loop integral of the subtracted GPDSE can 
only be recovered after taking into account the UV logarithmic behaviour for large-momenta dressing 
functions Q. Furthermore, /uv, behaving too as A 2 , should also be added in the r.h.s. of Eq. ((3T|) in 
order to write the renormalised GPDSE. Thus, the dependence on A but not the factor in front of it 
can be inferred from the GPDSE with only the information of the asymptotics for small-momentum 
dressing functions. 



2.1.5 The two classes of solutions 

The starting point for the following infrared analysis will be the Eq. (f3"Tj) for small A, where we will 
try to make the dependences on fc, k and A of the two sides match each other. 

The case ap ^ (scaling solution): We will first study the case ap ^ 0. Then, the l.h.s. of Eq. (|31l) 
can be expanded for small A as 

( A y~7 (41) 



F R (X 2 k 2 lf i 2 ) Fji(X 2 n 2 k 2 , fi 2 ) ~ I A(ii 2 ) 

and one obtains from Eq. (13T1) : 



N c gW) Z X AU 2 ) — ~ 1 , (42) 



where the dependences on fc, n and A of the numerator and the denominator should cancel against each 
other. Using for /ir the form given after Eq. (|3H|) . we find three possible situations: 



If otc + OiF > 1, applying Eq. (|4U)l in Eq. (|4"2"j) . we are led to the conclusion that only ap = — 1 (and 
ctQ > 2) satisfies this last equation and could be an IR solution for GPDSE. However, such a solution 
appears to be in a clearcut contradiction with the current lattice simulations. 
If cue + ap = 1, there is no possible solution because the logarithmic behaviour of Tir in Eq. (139p 
cannot be compensated by the powerlike one in the denominator of Eq. (1421) . 



and 



• If ac + dp < 1 , Eq. (|37|) combined with Eq. (|421) implies the familiar relation 2ap + uq = 
we have then: 

AT 2 ( 2 \ 7 (^ 2 )) 2 £(M 2 ) - / «G \ -i f,o, 

N c 9r{^ ) Z\ <py — ^-'"gJ - 1, (43) 

7 The multiplicatively renormalisable (MR) truncation scheme corresponds to letting /Iqcd — ► oo. Therefore, 
the scale qo being of the order of Aqcd, power laws with o>f + etc > 1 cannot be accepted as solutions of 
the GPDSE (se e, f or instance, [68|). The same argument holds also for of + ac = 1, because the ghost-loop 
integral in Eq. PTJ)) diverges as A — > for any go fixed as well as for qo — > oo for any fixed A. 
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An immediate consequence of this last condition is the freezing of the running coupling constant at 
small momentum. If the rcnormalization point, /k, is arbitrarily chosen to be very small in order that the 
dressing functions observe the power laws at k 2 = /i 2 , one obtains A(pi 2 ) = fi~ 2aF and B(/j, 2 ) = [i~ 2aG . 
Eq. (O then reads 

N c g 2 R (n 2 ) Zi (-|,a G ) ~ 16tt 2 , (44) 

and should be satisfied for any small value of /x. Consequently, it should remain exact as (i — > and 
provides the small-momentum limit of the running coupling (which is independent of the infrared 
constants for ghost and gluon dressing functions). 

In particular, if a G = 1, one has </>(— 1/2, 1) = 8/5 and thus 

N c g 2 R (p 2 ) Zi - IOt 2 , (45) 



The case ap = (decoupling solution): The case ap = is particular in that the leading contributions 
to the two occurrences of F in the l.h.s of cq. (f3"T|) cancel against each other. We have then to go 
one step further, taking into account the subleading terms. Defining Fir by means of Fj R (q 2 , /j 2 ) = 
A(n 2 ) + F IR (q 2 ,fi 2 ) we rewrite the l.h.s of eq. (23 as -(F IR (\ 2 k 2 , y 2 ) - F IR (X 2 K 2 k 2 , fi 2 )) / A 2 (fi 2 ) and 
use the known IR behaviour of Iir(A) from eq l3"U)0 in the r.h.s. of eq. (|3"Tj) to get 



(A^ 2 ) + A 2 (v 2 )q 2 lnq 2 ^ if a G = 1 
\ A([i 2 ) + A 2 {^ 2 ){q 2 ) a p ] otherwise 



(46) 



Furthermore, not only the subleading functional behaviour of the dressing function can be constrained 
but also the coefficient A 2 in Eq. (|4"o]) . In fact, if we plug this equation into the l.h.s. of eq. (|3"Tj) and 
expand we obtain : 



{A{^)f 



2 , „ ~ ffc 2 (l- K 2 )A 2 lnA 2 
No g R {n ) Z, / m (A) c | (A2fc2)Q ( 2 , (1 _ K2Q 



^2(^ 2 ) 

Let us consider now in more detail the three possible cases 
• If Q.Q < 1, we obtain from eqs. (|37l47p that 

(A{^ 2 )fB{^ 



if OLQ = 1 

(2) 

p ) otherwise 



(47) 



„( 2 ) 



a G 



Then, 



A2(v 2 ) 



NcgUn 1 ) Z l 0(0, a G ) 



where, according to eqs. (|34I36[) 0(0, oiq) is given by 



0(0, o G ) 



2a G (a G + l)(a G + 2)(l-a G ) 
Similarly if a G = 1, Eq. ([39]) applied to Eq. (J47|) leads to 

(A^ 2 )) 3 ^^ 2 ^ 



^2(M 2 



N C g R {^) Z 1 



64tt" 



(48) 



(49) 



(50) 



At last, if czq > 1, eqs. (gDJ) and (g7|) imply: 



(2) 



1 



a ghost dressing function which 



behaves quadratically for small momenta, In this case, however, as already said the ghost loop 
cannot be evaluated using the IR power laws over the whole integration range and it is therefore 
not possible to solve the GPDSE consistently, nor even to determine the small-momentum behaviour 
of the dressing functions, without matching appropriately those power laws to the UV perturbative 
formulas. Thus, we are not able to derive a constraint for the next-to-leading coefficient, A 2 {^ 2 )- 
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In summary, the GPDSE admits IR solutions with ap = and any ag > 0, provided that 



Fm(q 2 ,fi 2 ) 



A(^) l-0(O,a G ) 



16tt 2 



A 2 ^ 2 )B^ 2 ){q 2 y 



[A(^) + A 2 (n 2 )q 2 



■A 2 {n 2 )B(n 2 ) q 2 \nq- 



a G < 1 

OLG = 1 
a G > 1 



(51) 



where g 2 ^ 2 ) = iVc g 2 R (^ 2 )Zi and 0(0, a G ) is given in Eq. ([49]) . The gluon dressing function is supposed 
to behave as indicated in Eq. (|27|) . In particular for ag = 1, the gluon propagator takes a finite 
(and non-zero) value at zero momentum, B(fj, 2 ), after applying MOM renormalisation prescription at 



2.2 The low-momentum solutions and the gluon mass 

As will be seen in the next section, the current lattice data strongly supports a decoupling solution 
which does not obey 2ap+ac = and in which a G = 1- Furthermore, lattice data can also be very well 
accommodated within DS coupled equations in the PT-BFM scheme [ltl [2(| and within the so-called 
refined Gribov-Zwanziger approach [42j, leading in both cases to decoupling solutions for gluon and 
ghost propagators. Then, as Eq. (|51[) reads for, this implies a ghost dressing function proportional to 
q 2 logq 2 while the gluon propagator takes a constant at vanishing momentum (ag = 1) or, in other 
words, is "massive" . 

Indeed, it is well known that the Schwinger mechanism of mass generation (69| can be incorporated 
ino the gluon propagator DSE through the fully-dressed non-perturbative three-gluon vertex and gives 
rise to the generation of a dynamical gluon mass such that (l9l : l70l | , 

D-\q 2 )~q 2 + M 2 (q 2 ) . (52) 

In particular, it has been shown that a power-law running mass, 




^v) = =Fh* : :j ' > (53) 



appears as a solution for the coupled ghost and gluon propagator DSE in the PT-BFM truncation 
scheme (7l| and, also, as a consequence of the dimension-4 gluo n condensate in the OPE expansion of 
the gluon self-energy in the Pinching Technique framework [72j . 

Having this in mind, the authors of ref. [35| applied the following simple model 0: 

2 2> _ Gie(« 2 ,^ 2 ) B(^) B(/. 2 )/, 5 2 ,„(q' 



for a massive gluon propagator, in order to compute the O(g 2 )-correction for the low-momentum ghost 
dressing function and then prove that this low-momentum behaviour is controlled by that gluon mass 
and by the zero-momentum value of the effective charge defined from the Taylor-scheme ghost-gluon 
vertex in ref. j29j . 

The work of ref. (3|| can be easily overviewed if we re- write Eq. (|2~4"|) as follows 
1 1 9rNc f ,4 F R (q) f 1 ^ (k-q) 2 -k 2 q 2 

= N C9 2 R (I m (k 2 ) - I m (p 2 )) (55) 

8 This is a renormalized massive gluon propagator, as given by Eq. (|52|l . where the gluon running mass 
appears to be approximated by M(q 2 ) ~ M(0) = mo = M. This is, for instance, a very good low-momentum 
app roxi mation for the running mass given by Eq. (|53|) with mo ~ 0.5 GeV and Aqcu ~ 0.3 GeV (see fig. 11 of 
ref. HI). 
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for the two momenta fc and p, which will be taken to be small (as compared with M and qo, the 
cut-off separating IR and UV modes in the ghost-self-energy integral), and where the ghost-gluon 
transverse form factor, Hi, has been explicitly replaced by 1 and the dependence on the renormalization 
momentum, /i 2 , omitted by simplicity. Then, if one inserts the form Eq. (|54[) of the gluon propagator 
into the integrand and focuses on the low-momentum behaviour of the ghost dressing function for the 
decoupling case (of = 0), the ghost dressing function will be replaced by the constant leading term, 
A(/j, 2 ) from Eq. (f2"T)l . and one will obtain: 

i i 



k 2 + q 2 - 2kq cos(6>) M 2 + k 2 + q 2 - 2kq cos(0) 

1 A(fi 2 ) B(/i 2 ) f qo ^ / ., s 3k 2 -q 2 ,,3g 2 
q 3 F(q) dq 6(k - q) — ±- + 6(q - k)- 



32 7T 2 M 2 J y w ^ V ki I 4 

6k 2 q 2 (M 2 + fc 2 + q 2 ) - (M 2 + fc 2 + q 2 f + ((AJ 2 + k 2 + g 2 ) 2 - 4fc 2 g 2 ) 3 / 2 

2fcV 



(56) 



The integral given by Eq. (|56p can be first analytically obtained and then expanded up to the first 
order in k 2 /M 2 to give: 

2 1 A( M 2 ) j?(/x 2 ) 
/lR(fc) = _ 32^ M 2 



3 „ r2l ,^ 2 + 90x fc \ , 11,2 fc2 , /, M \ k 2 q 2 M 2 

- 2 M*lo g{ ^^) + - log( — ) - -fc 2 + - log(l + ^) - y 



from where one needs only to keep 



„9n 1 A(u 2 ) Bin 2 ) , , fc 2 . 11, 9 



64 7T 2 M 2 V by M 2 ' 6 

because the terms in 0(M 2 /q 1 ) are neglected, while the constant terms and the ones logarithmically 
divergent (~ log (qo)) happen to be cancelled when applying the subtraction for two different momenta 
in the r.h.s. of Eq. (1551) . 

Thus, if one replaces Iir by the result of Eq. (j58|) in the r.h.s. of Eq. ((55)) . the subleading term for 
the ghost dressing function should be: 



2 



1 + w St(0) w 



\ q 2 11 



where: 



5r (0) = Jim (q 2 + M 2 ) = M 2 9 -^F 2 R (0, M 2 ) D lK (0, » 2 ) 

= S -M1a 2 (h 2 )B(^ 2 ) . (60) 

Here olt = g^/(ATr) is the perturbative strong coupling defined in the Taylor scheme [zH, while ax is 
the non-perturbative Taylor effective charge defined as an extension of the Taylor ghost-gluon coupling 
in ref. (29|, similarly to the way the "pinching technique" (PT) effective charge was from the gluon 
propagator in ref. |74l |. 

In the forthcoming sections [3] and [H Eq. ([5T)]) will be confronted to lattice and numerical DSE 
estimates for the ghost dressing function in the low-momentum domain. 



9 The "pinching technique" [7(| implies a resummation of the diagrams for the perturbative gluon propagator 
expansion leading to a redefinition of the propagator such that it observes a QED-like Ward identity, thus 
providing us with a way to construct an IR effective charge as happens in QED. 
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2.3 Some constraints from Slavnov- Taylor identities 

In the previous section, we have analysed the infrared behaviour of the GPDSE solutions and found 
that the ghost dressing function can cither diverge at vanishing momentum (otp = — ciq/2 with uq > 0) 
or take a finite value (op = with any olq > 0). As appendix iBl shows, the GPDSE themselves can 
be derived from the general Ward-Slavnov- Taylor equation [75| . The Ward-Slavnov- Taylor identities 
(WSTI) can be derived formally from the gauge invariance of the path integral, eq. (|118[) of appendix 
IBl as shown in [76j ], This is the case in lattice simulations. If the path integral is limited to a domain 
of the configuration space such as in the Gribov-Zwanziger approach, used in [42| . the STI may not 
be satisfied. We assume a gauge invariant path integral and will now invoke the WSTI for general 
covariant gauges relating the 3-gluon, F\ pv (j>^ q, r), and ghost-gluon vertices, 



F{P 2 ) 
G(r 2 ) 

F(p 2 ) 



P\r\^(p,q,r)-——(5 p „r 2 - r p r„)r pp (r,p; q) 



(61) 



G{q 2 ) 

to shed some light on that matter [77l [78| . Using for the ghost-gluon vertex the general decomposi- 
tion [ZSl 

r^, p (p,q;r) = 5 vfl a(p,q;r) - r v q p b(p,q;r) + p v r p c(p,q;r) 

+ r v p p d(p,q;r) + pvp^ e(p,q;r) , (62) 

and multiplying by r„ both sides of Eq. (|61[) . one obtains: 

Ftp 2 ) 

r v pxr Xpu {p 7 q,r) = X{q,p;r) [(q ■ r)g M - q 2 r p ] ; (63) 

where 

X(q,p;r) = a(q,p;r) - (r ■ p) b(q,p;r) + (r ■ q) d(q,p;r) . (64) 



Since the vertex function, _T, in the l.h.s. of Eq. (|63[) is antisymmetric under p r and A <-> v, one 
can conclude thatfl: [77^: 

F(p 2 )X(q,p; r) = F(r 2 )X(q, r;p) . (65) 

This last result is a compatibility condition required for the WSTI to be satisfied that does not involve 
the 3-gluon vertex and implies a strong correlation between the infrared behaviours of the ghost-gluon 
vertex and the ghost propagator. Now, under the only additional hypothesis that those scalars of 
the ghost-gluon vertex decomposition in Eq. (|62p which contribute to the scalar function X defined in 
Eq. (|64p are rcgulai0 when one of their arguments goes to zero while the others are kept no n- vanishing, 
one can consider the small p limit in Eq. (|65|) and obtain: 



F(p 2 )X(q, 0; -q) = F(q 2 )X(q, -q; 0) + 0(p 2 ) (66) 

This has to be true for any value of q, which implies that F(p 2 ) goes to some finite and non-zero 
value when p goes to zero, since neither X(q,0; ~q) nor X(q, — q;0) are presumably zero for all 
values of q. Rephrased in terms of infrared exponents, the latter argument implies that cxf = 0. 

To reach the above conclusions we did not appeal to the properties of the 3-gluon vertex, apart 
from the symmetry under the exchange of gluon legs. If one assumes in addition that the longitudinal 
part of the 3-gluon vertex also behaves regularly when anyone of its arguments goes to 0, the others 

10 We work, of course, on the energy-momentum shell, so that the relation p + q + r = holds 

11 Note also that, for our purposes, it will actually be enough to restrict, but not forbid, the possible presence 
of singularities in the scalar coefficient functions provided that they could be compensated by kinematical 
zeroes stemming from the tensors. 
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being kept non- vanishing, a divergent gluon propagator at vanishing momentum will be obtained from 
Eq. ^2|s3|23: 



lim p\r u r Xflu {p,q,r) 

q— >0 



rxr v r\ MV (-r,Q,r) 

«2 



lim 

q->0 



Hp 1 



G{q 2 



implying that 



(q-r) 



% X(q,p; 



(67) 



%W 2 



= 



(68) 



Of course, as far as it involves a vertex with longitudinal gluons which have not been very extensively 
studied, this last conclusion is not as clean as the previous one about the ghost dressing (according to 
authors of rcf. (8lj a soft kinematical singularity appears for the landau-gauge 3-gluon vertex, however it 
does not concern our proof relying on the regularity of the longitudinal-longitudinal-transverse 3-gluon 
vertex) . 

In ref. [77j . we showed that only a very mild divergence, for example of logarithmic type, could be 
compatible 

(although very unlikely) with current LQCD results for the gluon propagator. The IR analysis of 
the previous section can be straightforwardly extended to this case by generalizing 



G m {q\^) = B(u. 2 )(q 2 )^ \og»[- 



(69) 



the effect of which is to modify Eq. (|51l) with 



iW,M ) 



A(^) (l - 4>(0,a G )^^A 2 (^)B(^)( q 2 r° log- (<T 2 )) a G < 1 

A(^) (l - { X\)L* A V W) I 2 W" +1) (I- 2 )) <*G = 1 ( 7 °) 
{A(^)+A 2 (^)qHog v (q- 2 ) 



a G > 1 



where only the power of the logarithm is then modified. 

Sticking now to the case where cxf is zero (for the reasons explained above) and a G is 1 (as suggested 
by the lattice results) we are left with 



F m {q 2 ,^ 



1 



(v + 1)64tt 2 



F m (p,nyB(ji*) q 2 log 



rl 2 



(71) 



according to whether there are logarithmic corrections to the gluon propagator (v ^ 0) or not {v = 0). 
Here, M is some scale which is out of the scope of the IR analysis we performed in the previous section 
and, if v = 0, -B(/i 2 ) = G|^(0, is the gluon propagator at zero momentum. 



3 Low-momentum Green functions lattice results 

As was recalled in the previous section, the mechanisms usually invoked to explain the confinement 
imply specific behaviours in the infrared for the Green's functions of the theory : 

• A sufficient condition for the Kugo-Ojima criterion^ to be satisfied would be the divergence of the 
ghost dressing function, F^ 1 (q 2 ) ~ as q 2 -> (see ref. [60|; l6l| ) 

• The Gribov-Zwanziger scenario implies the vanishing of the gluon propagator, i.e. G(q 2 )/q 2 ~ as 
q 2 -> (ref [12) 



Th e Kugo-Ojima scheme also implies that there be no massless pole in the transverse gluon propagator 
(cf [E0]), a condition which is weaker than the vanishing advocated in the Gribov-Zwanziger scenario. 
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In addition, the dressing functions should obey Dyson-Schwinger equations, from which it was 
inferred that the ghost and gluon infrared exponents (cf. equation ^ above) should satisfy the relation 
aa + 2a.F =0 (referred to in the following as Rela) Accordingly Zwanziger [59| and Lerche and von 
Smekal [66| predicted a value olq = —2ap ~ 1.19 

Since those results rely on theoretical conjectures and since using the Dyson-Schwinger equations 
demands an unescapable and not fully controlled truncation of their infinite tower, the technique of 
lattice simulations has become an alternative and totally independent means to get reliable model- 
independent information on the Green's functions. The first attempts to measure propagators took 
place at the end of the eighties for the gluon [H; [13 and in the second half of the nineties 
for the ghost [HI; H3; HI)]. Since then, thanks to the huge increase in the performance level of the 
computers, it has become possible to reach larger volumes while keeping a small enough lattice spacing. 
This circumstance is essential to get an ever more detailed insight in the infrared behaviour of the 
propagators and vertices. 

The simulations have been performed using a variety of setups : 

• choice of the gauge group (SU(2) or SU(3)) Theoretical arguments lead to the conclusion that the 
qualitative features of the IR Green's functions should be independent of this choice [HI; H(| . 

• dimensionality d of the space, ranging from 2 to 4. The relation (Rela) is actually the restriction 
to the d — 4 case of a more general e?-dcpendent one. 

• choice of the gauge action : either standard Wilson or improved versions, quenched or unquenched. 

• lattice geometry (isotropic or not), spacing, lattice size. 

We present in this section an overview of the results of these simulations, with emphasis on the 
case of an SU(3) gauge group and of the pure Yang-Mills theory (quenched QCD). The case of the 
SU(2) gauge group has been considered very thoroughly by Cucchieri and Mendes (see for instance 
[89j and the references therein). A comparison of the dressing functions for SU(2) and SU(3) has been 
performed in ref. (25[ and shows that the dressing functions are remarkably close for both the gluon 
and the ghost, in accordance with the theoretical expectations. 



3.1 Ghost and gluon propagator results in Landau gauge 

We summarize in table [T] the results concerning the IR properties of the propagators which have been 
reported in the literature. They are usually, but not always, given in terms of values of the exponents 
and «g- The reader should keep in mind that in many of the latter cases, the authors have attempted 
to describe both the gluon and the ghost propagator with only one parameter, automatically taking 
into account the relation (Rela), but simultaneously introducing some kind of bias in the fit, all the 
more as, as we shall see, this relation is not satisfied by the data. Actually such a one-parameter fit 
revealed impossible in several cases [9(J ■ 

While the first measurements on the lattice resulted in a value of ag compatible with 2, i.e. with 
a gluon dressing function behaving as k A at small fc's [Ifl, the table shows that there is nowadays a 
general agreement (see also ref. [89( ) in favour of a solution incorporating : 

1. a gluon dressing function going to like q 2 , leading to a propagator remaining finite and non-zero 
( "massive gluon" ) 

2. a ghost dressing function going to a non-zero constant, i.e a ghost propagator behaving as 1/q 2 
("free ghost") 

Recent numerical data (regarding the quenched case ) on very large lattices are to be found in |105| . 
They are visualised in figurc[T]bclow, in which the left panel is borrowed from [l05l |. Similar simulations 
are now underway for the unquenched cas. The preliminary results appear to be qualitatively compati- 
ble with the general picture we have described. We defer the discussion of the different artefacts which 
can affect the results to a special subsection but, meanwhile we can only make ours the statement of 
ref. : " The current paradigm is that of a massive gluon and a free ghost" . 



19 



Ref 


year 





Lat. Size 


gluon 


ghost 


[91; 92] 


2001 


5.7 


16 3 x 32 






-id- 


-id 


4.38 (+) 


16 3 x 32 






-id- 


-id 


3.92 (+) 


10 3 x 20 


finite 




-id- 


-id 


3.75<+> 


8 3 x 16 


non zero 


n.m. 


-id- 


-id 


3.92 (+) 


16 3 x 32 


D(0) 




-id- 


-id 


4.1<+> 


12 3 x 24 






-id- 


-id 


6.0 


32 3 x 64 






[93] 


2005 


5.8 


16* - 32 4 


D(<f) 




-id- 


-id 


6.0 


16 4 ,24 4 , 


decreases 


a F ~ .25 


-id- 


-id 


6.2 


16 4 ,24 4 


with q 2 




[24; 94J 


2005 


5.75 


32 4 


a G = .864(16) 


a F = -.153(22) 


[95; 96] 


2006 


6.0 


(8 a -,18 3 ) x 256 


a G = .996 - -1.05 


n.m. 


[97] 


2006 


5.8 


24 4 ,32 4 


finite 


a F ~ .2 


-id- 


-id 


6.0 


16 4 - 48 4 


non-zero. 


or 


-id- 


-id 


6.2 


16 4 ,24 4 


D(0) 


log- like 


[22] 


2007 


5.7 


56 4 - 80 4 


-id- 


not power- like 


[98] 


2007 


6.0 


(8 3 - 18 3 ) x 256, 






-id- 


-id 


-id 


16 3 x 128,48 4 ,64 4 


a G = 1.07(28) 


not power- like 


-id- 


-id 


6.2 


64 4 






[99] 


2008 


6.0 


16 4 - 64 4 , 


not 




-id- 


-id 


-id 


(8 3 - 18 3 ) x 256 


conclusive 


n.m. 


[100; W1J 


2009 


6.0 


16 4 - 80 4 


not 


n.m. 


-id- 


-id 


5.7 


8 4 - 44 4 


conclusive 




[102: 103: 1041 


2009 


6.0 


32 4 






" -id- " 


2009 


5.8 


20 3 x 32 


non zero 


n.m. 


-id 


2009 


6.0 


16 3 x 32 


D(0) 




[105] 


2009 


5.7 


64 4 - 96 4 


a G = 1. 


of = 0. 



Table 1 Summary of the infrared behaviour of the gluon and ghost propagators from lattice simulations, 
restricted to the SU(3) case; the infrared exponents, a G and a F , are given when available; "n.m." stands for 
"not measured". The "( + )" mark refers to the use of the Liischer-Weisz improved action (cf. infra). 




Fig. 1 Gluon propagator (left, borrowed from |105lj ) and ghost dressing function (right, from [621] ') 



3.1.1 Lattices vs Dyson-Schwinger equations 

The fact that the most recent results of lattice simulations, ac = 1 and ac = 0, were actually 
compatible with the Dyson-Schwinger equations was first demonstrated in reference [26^ by a numerical 
analysis of the GPDSE where the gluon propagator lattice data were plugged into the kernel and the 
transverse form factor was supposed to be constant for all momenta. We will come back to this result 
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in sec. HJ where the agreement of lattice data and the numerical DSE results for the decoupling case 
will be manifest from fig. [7] 



P=5.7, N=80 
P=5.7, N=64 
|3=5.8, N=32 
P=6.0, N=24 
M=0.50 GeV 




LO 



k (GeV) 



Fig. 2 The gh ost dressing function from lattice data [621 ; I22l| is pretty well described by the low-momentum 
formula of Eq. (|59|) with M = 0.50(2) GeV, as shown in ref. [351] from where we borrow the figure. The curve 
for the low-momentum formula is drawn as solid black line inside the fitting region and as dashed one outside. 



Furthermore, as we discussed in the previous section, in the subsection !2.2[ the authors of ref. |35j 
provided us with a low- momentum analytic expression for the ghost dressing function, Eq. (1591) , derived 
from the asymptotic analysis of the GPDSE. This result was also succesfully confronted to the lattice 
data for the ghost propagator dressing, in particular with those obtained from very large lattice 

from which we borrow Fig. [5] where 



simulations in ref. [22]; this comparison was performed in 
the ghost propagator lattice data are shown to behave pretty well as Eq. (|55|) asks for with a gluon 
mass, M = 0.50(2) GeV, in the right ballpark (roughly from 400 MeV to 700 MeV) define d by 
phe nomenological tests [l06j ] or direct lattice measurements from the gluon propagator [9l|; Il02t 



3.2 The ghost-gluon vertex 

Among the numerous possibilities to define the QCD renormalised coupling constant gn a particularly 
attractive one, for use in lattice simulations, is based on the ghost-ghost-gluon 3-point function together 
with a specific MOiVf-type renormalisation scheme in which the incoming ghost has zero momentum. 
We have intoduced in eq. (|T8|) the 2 form-factors Hi and H 2 of the vertex; in terms of those factors 
the generic definition of gn would be 

g 2 R {H 2 )= lim g 2 (A 2 )Z 3 ^ 2 ,A 2 )Z 2 ^ 2 ,A 2 )(H 1 (q,k;A) + H 2 (q,k;A))\ q 2 = ^ (72) 

A— >oo 

In the specific kinematical situation we have just mentioned, Taylor's non renormalisation theorem 
states that Hi(q, 0; A) + H 2 (q, 0; A) = Zi(p 2 ) = 1 from which follows that 

a T {^ 2 ) ee 9 ^fl = lim 9 M1 G ^ 2 1 A 2 )F 2 ^ 2 1 A 2 ). (73) 

47T A-^oo 47T 

These lattice data were also obtained with the SA gauge-fixing algorithm so as to deal a s w ell as possible 
with the Gribov ambiguity but, for instance in determining the gluon mass by fitting Eq. (|59|) to the ghost 
dressing data, some systematic uncertainty should be admitted to come from the gauge-fixing procedure. 
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The remarkable feature of Eq. (|73l) is that it involves only G and F so that no measure of the 
ghost-gluon vertex is needed for the determination of the coupling constant. 

It is thus particularly suitable for use on the lattice, since measuring it does not demand any delicate 
3-point fu nctio n computation ; therefore it has been extensively advocated and studied (see for instance 
reference |108| ). Were the relation (Rela) fulfilled, this quantity should go to a finite non-zero limit in 
the infrared. What is observed on the lattice is quite different from this expectation : the data invariably 
show a fast decrease to 0, as displayed in figure[3] (see also Fig.|S]in the next section H]). The mechanism 
which leads to the apparition of those 2 types of solutions has been analytically discussed at length in 
the second section. Here we just want to stress the fact that the lattice simulations undoubtedly favor 
the "decoupling" case. Of course this argument does not say anything on aj? and ag separately but 
it shows, at least, that one should use 2 independent exponents to describe the infrared behaviour of 
the propagators. 

To our knowledge, the direct measurement of the ghost-gluon vertex that has been performed with 
the SU (3) gauge group appears to be very noisy and helps to conclude n othin g beyond the fact that 
the transverse form factor, Hi, in Eq. (|72p is pretty close to 1, as refs. j90t Il09{ clearly showed for the 
case of a ghost-glu on v e rtex with vanishing incoming gluon momentum. There exist also simulations 
in the SU{2) case [llOt llllj . which constitute a very clear direct check of Taylor's theorem. On the 
other hand, it is not clear whether the data presented in this reference for a s favour the scaling or the 
decoupling solution. 
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Fig. 3 Recent and very accurate data for ar(q 2 ) borrowed from ref. [1051] (figure 5), which shows a low- 
momentum behaviour clearly compatible only with the decoupling solution. 



Furthermore, as will be discussed at the end of app. [K\ Eq. ([73")) missed some non-perturbative 
corrections which, although playing no role for the running in the UV domain, need to be taken into 
account when integrating numerically the GPDSE. 



3.3 The 3-gluon vertex 

It is of course also possible to define the stron g coupling constant directly from the three gluon vertex. 
Let us recall shortly this definition of a s (p 2 ) jl!2t Ill3| . We consider the three-gluon Green function 



G(3) uiU3Us(Pl5 2 J 2,Z>3) at the symmetric point, p\ = p\ = p\ = u 2 
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The tree-diagram three-gluon vertex is given by g s T tree with T tree defined by 



T, 



[Vi/*i(Pi -Pi)n' 3 +cycl. perm.] ] [ 1(5 



i=l,3 



(74) 



The three-gluon Green function may be expanded on a basis of tensors. We are interested in the scalar 
function G^ 3 ' (fi 2 , fi 2 , fi 2 ) which multiplies T tree . It is obtained by the following contraction 



_^ jaia 2 



(Pl,P2,P3) 



18^ 2 24 

(Pi ~P2) M3 (P2 ~P3)mi(P3 -Pl) M 2 

The Euclidean two point Green function in momentum space writes in the Landau Gauge: 



where ai,a2 are the color indices ranging from 1 to 8. 

Then the renormalised coupling constant is given by |112| 



9r{^ ) 



g^HpIpIpDz 3 / 2 ^ 2 ) 



G(2)(p2) G (2) ( p2 )G (2) (p 2 



where 



Z 3 ( M 2 ) = M 2 G< 2 V). 



(75) 



(76) 



(77) 



(78) 



This method has been used in the infrared regime in ref. [1141 ] . It leads to a coupling constant which 
can be nicely fitted by a p 4 law compatible with an interpretation in terms of an instanton gas, as is 
shown in figure 2] 



3.4 Reflexion positivity violation 

Although the link between the confinement property and the infrared behaviour of the propagators is 
not fully understood, it is generally admitted that the spectral functions of the latter should not be 
positive definite, since, in the coloured sector, no physical posit ive norm state can contribute. That 
this is actually the case has been verified by several groups [13; ll 1 5t Il03l ; [l04; 102]. To perform this 
check one considers the quantity : 

fJ,,X 

The 0- momentum gluon propagator can be obtained from C by integrating it over t. Therefore C(t) 
must change sign in order that D(0) vanish, but this necessary condition is of course not sufficient. 
Note that C(t) is similar to what one considers in hadronic physics simulations to "measure" the 
masses (typically the pion mass). It had already been noticed by Bernard et al. ([Ifl) that the effective 
mass m(t) = Log(C(t + l)/C(t)) increases with t, in contradiction with what happens in the "p hysical" 
situation where only the lightest state survives as time increases. Sternbeck an his collaborators 97l : lll5| 



have sh own by d i rect inspection that C(t) becomes negative for large enough As for the japanese 
group ( [lQ2t Il03l ; Il04| ). performing a fit to the lattice data and computing the spectral density o f the 
fitting function they show that this density is almost everywhere negative. Cucchieri and Mendes jl 16| 
consider the gluon propagator in x-space as a function of the 4-dimensional distance and also observe 
that it becomes negative as the distance increases. Thus the four methods reach the same conclusion, 
namely that the gluon propagator as measured in lattice simulations does violate the reflexion positivity. 

14 They have furthermore checked that this remains true in the unquenched case. 
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10000 




Fig. 4 Coupling constant in a cooled lattice (L = 24, /3 = 6.0) after 200 cooling sweeps. The solid line 
corresponds to the fit discussed in the text. The horizontal axis is given in GeV, assuming for simplicity the 
lattice spacing of the thermalised configurations, a -1 = 1.97 GeV. 



3.5 The artefacts 

3.5.1 Finite spacing effects and rotational invariance violation 

On the lattice, in dimension 4, the usual 0(4) symmetry which prevails in continuum euclidean field 
theory is broken down to the hypcrcubic H{4) invariance. As a consequence the number of independent 
scalar invariants is increased to 4 : = X}„p^, 71 = 2, 4, 6, 8. Evidently, dimensional arguments show 

that the higher invariants should actually appear as a 2 pM , a 4 p^ • • • so that they are also higher order 
corrections in terms of lattice spacing. This manifests itself through the "fanning" or "fishbone" effect : 
when displaying a scalar function versus p one observes a series of characteristic fringes due to the 
dependence on the higher invariants. It makes the signal noisy at large momenta and one has to cure 
this effect in order to identify the "physical" value. The simplest way to realize this goal consists in 
keeping only the momenta which minimize p^ as compared with pl 2 l; in practice this amounts to 
keeping only those momenta the components of which differ at most by a, few units from the diagonal 
case (cf [H(). A more efficient technics has been devised later on (l 1 7t 1 1 18t 1 1 19| . The idea is to perform 
an expansion of the scalar functions into a series ot the invariants : 

fipWpWpWpW) = /(y 2] ,0,0,0) +pW^(pP] >0 ,0 > 0) +•■■ (79) 

The first term of the series is the desired continuum scalar function, up to rotationally invariant 
finite lattice spacing effects. The latter are expected not to be important in the deep infrared domain 
since in this region pa -C 1. At moderate momenta the situation regarding those discretisation effects 
depends on the quantity under scrutiny ; for a given physical momentum a comparison between data 
at /3 = 5.8 and f3 = 6.0 shows a perceptible but moderate increase of the gluon dressing function G 
while the ghost dressing function seems unsensitive to the effect (93|. 

3.5.2 Finite volume effects 

The size of the lattice is, for several reasons a very important parameter. 
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First of all it determines how deep in the infrared it is possible to go. The 0-momentum value of 
the gluon propagator and of its dressing function can be obtained directly on the lattice. Of course the 
zero-momentum value cannot be used directly to fit aa, since including it in the data would force it 
to be equal to 1. But it is useful in making it possible to compare the data to possible analytic forms, 
or to look for possible discontinuities, provided sufficiently small momenta are available. The situation 
is different for the ghost: in this case the zero momentum is not accessible. Then the determination of 
the ghost infrared exponent relies crucially on the possibility of getting as close as possible to ze ro. 

Second, the DSE equations on a torus in the continuum have been studied in reference [120| |. The 
authors compare numerical solutions for various volumes to the infinite volume one and show that a 
minimum volume of (10 — 15/m) 4 is necessary in order to observe the onset of the infrared behaviour 
of the Green's function and that only much larger ones could allow a reliable estimate of the infrared 
exponents. 

A third reason why volume effects might be important has to do with gauge fixing and the prob- 
lem of Gribov copies. It is believed that, as the volume increases, the integration measure over field 
configurations should concentrate on the boundary of the Gribov region. This entails that the ghost 
propagator must increase with the volume, because the smallest eigenvalue of the Faddeev-Popov 
operator is located on the boundary and goes to in the infinite volume limit. 

A number of strategies have been elaborated in order to determine whether or not the measured 
lattice quantities approach the physical values. The most direct one consists evidently in measuring 
-D(O), agandaf on a set of lattices of increasing volumes, to study their volume dependence and, if 
possible, to determine their infinite volume limit, either first determining for each finite volume the 
infrared exponent, then extrapolate to get otoc or extrapolating the raw data for the propagators and 
then determine direct ly oioo . Others use a set of V-dependent bounds to constrain the IR behaviour of 
the propagators [l2lj . Finally, in ref. (98| . the authors have developed a method which minimizes the 
volume effects while checking the power-like dependence. 

The situation is not completely settled yet, but : 

1. The sizes mentioned above as minimal requirements to obse rve the switching to the infinite volume 
regime have now been passed. Volumes of (16/m) 4 in [l05l | and up to (27/to) 4 for SU(2) in [89| ) 
have bee reached wit hout any sign of a change in the curves showing up. 

2. While some authors [1051 : 192^ present results according to which the small momentum gluon prop- 
agator is only slightly dependent o n the volume and conclude to a finite infinite- volume limit, 
according to others (see for instance [lOlj ) the present data are compatible with both a finite and 
a vanishing popagator. 

3. A similar situati on pr evails for the ghost propagator : Bogolubsky et al. [l05l | conclude to finite 
F(0), but others [lOll ] consider that F cannot be described by a simple power law. 

Altogether the curves for G show a remarkable consistency between the small momentum plateau 
and the direct 0-momentum measurement. As volumes are increased, less and less room is left for a 
possible bending down. Then, taking also into account the very neat results of the SU(2) study of 
Cucchieri and Mcndcs on very large lattices we are led to decide in favour of the decoupling solution, 
although more work is still needed to be absolutely sure that neither G(0) ^ nor a finite F(0) are 
finite volume artefacts. 

3.5.3 Lattice anisotropy effects 

A means to reach small values of the momentum without having to deal with prohibitively large lattices 
consists in using anisotropic lattices with a very large number of nodes in one of the directions (usually 
the temporal one) . The question whethe r thi s asymmetry might induce artefacts has been raised in the 
SU(2) case by Cucchieri and Mendes ( |122| ) who concluded that, although the IR behaviours of the 
Green's functions looked qualitatively similar, they were quantitatively dependent on the geom etry ; 
later on it has been extended to SU(3) and considered in detail by Oliveira and Silva ( jl23l |). who 
compare, in the anisotropic case, the purely temporal momenta to spatial ones, as well as the isotropic 
and anisotropic situations for various types of cuts. Neither of these comparisons show any effect 
in the infrared region. However the same authors, but using a different approach (Cucchieri-Mendes 
bounds) conclude in [99^ that the data from symmetric and asymmetric lattices are incompatible, 
acknowledging at the same time that there is no theoretical explanation for the phenomenon. The 
situation still demands clarification. 
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3.5.4 Influence of the action 

To conclude this section we mention that the possibility that the choice of the action might introduce 
artefacts has been considered by Bonnet et al.[92[. They compare the results obtained for the gluon 
propagator from the standard Wilson action with the outcome of the Liischer-Weisz 0(et 2 )-improved 
action in nearly similar conditions of lattice spacing and volume. The raw data in the second case 
exhibit considerably reduced finite spacing artefacts, as was to be expected. Once the data have been 
treated for those a rtefa cts (see above in subsection 13.5. ip the agreement is excellent. This study has 
been extended in [124| to the unquenched case. Again, the gluon propagators obtained under the 2 
hypotheses present the same qualitative infrared features. 

3.6 The Gribov problem 

As was discussed in sec. ll.3[ Gribov (3(| first realised that in a non-Abelian gauge theory there remains 
a gauge ambiguit y ev en if imposing, in the case of Landau gauge, the constraint d^A a ^ = 0. This is also 
true on a lattice jl25{ where it is referred to as the "lattice Gribov ambituity" . On the lattice Landau 
gauge is fixed by minimising the functional 

F u [g} = -J2U[Tr(U^x))}= f o% g\ c? A\x) + 0{a A ). (80) 

This verifies d^A'^ = since d^A^ is the derivative of the functional in Eq. (jSOJ) as a function of 
the infinitesimal gauge transformations. Picking a minimum implies that we choose a gauge such 
that the Faddeev-Popov operator is positive (which is Gribov's prescription) since the Faddeev-Popov 
operator is the second derivative of the functional in Eq. (|80[) as a function of the infinitesimal gauge 
transformations. But even this restriction, which eliminates local maxima and saddle points, is not 
sufficient because there are many local minima. It was suggested to choose the absolute minimum on 
the gauge orbit [37}, however this is numerically out of reach. A reachable method was used [23|; [8?} 
consisting in taking a number of random copies, selecting the "best copy" which provides the minimum 
of the functional in this sample and compare it to the "first copy" appearing in the random sample. This 
allows to estimate the effect of decreasing the functional and the amount of scattering of the considered 
quantity. It has been applied to the gluon and ghost propagators. The general conclusion is that there 
is no dependence on the Gribov copy except for the smallest momenta (infrared). In that region the 
gluon propagator depends moderately on the copy, while the ghost does depend significantly : the 
ghost propagator is smaller for the best cop y tha t for the first copy. 

Moreover, a further SU(2) investigation |105| based on the application of the so-called "simulated 
annealing algorithm" (SA) gauge-fixing algorithm leads to the same conclusion: lattice results for ghost 
and gluon propagator behave as expected for a decoupling solution. This SA is a "stochastic optimiza- 
tion method" allowing quasi-equilibrium tunneling through functional barriers, with a stati stica l weight 
which is ex exp — F(j[g]/T where T is a "temperature" that should be taken to decrease [l26j . Then, 
in principle, with the appropriate T decrease and number of cooling steps, the SA requires the gauge 
fixing functional to take the extrema as arbitrarily close to the global extremum as wished; in other 
words, SA permits to select the "best copy" by fixing the Landau gauge in the "fundamental modular 
region" . Although it has been recently guessed that finding the Gribov copies as close as possible to the 
global extremum may maximally enhance t he in frared asymptotics of the ghost dressing function, as it 
would correspond with the scaling solution jl27l | , the comparison of S A results with the ones obtained 
by applying the standard g auge - fixin g method based on the "over-relaxa tion" algorithm seems not to 
support such a conjecture [1281 : Il29j . In particular, the authors of ref. |l29l | found the Gribov copy 
effect not to have any impact for momenta above a given p m j n depending on the physical lattice size 
(Pmin decreases when lattice size, aL, increases) and concluded that the SU(2) gluon propagator was 
compatible with a decoupling solution. These authors compared the gluon propagator results at any 
fixed momentum obtained by fixing the Landau gauge by the random choice of a first Gribov copy 
(fc) with the ones obtained by the choice of the best copy (be) as the extrema of the gauge fixing 
functional. They concluded that the discrepancy clear ly tends to disappear as the lattice volume in 
physical units increases (see the figs. 2, 3 and 4 of |l29l |). Thus, although some gauge-fixing ambiguity 
for the very low-momentum gluon propagator data is reported, it is claimed either to disappear in the 
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large-volume limit or to be avoided in practice by applying the SA algorithm to fix the Landau gauge 
and the best Gribov copy. 

On the other hand, the continuum limit of one Gribov copy is impossible to perform, but some 
statistical quantities re lated to these copies can be studied. This has been performed on the SU (2) non- 
Abclian gauge theory [l30l | . It was found that the probability to find a second copy strongly depends 
on the size of the lattice in physical units. There is a fast transition between high probability and 
small probability at a lattice length of about 2.75/ \fa where a is the asymptotic string tension and 
is estimated in QCD to be y/a ~ 0.5 GeV. This implies a critical length above which Gribov copies 
become frequent of ~ 1 fm. It is not surprising that in the small volume limit, close to the perturbative 
regime, Gribov copies tend to disappear. They are typically a non-perturbative phenomenon. That their 
existence depends on the physical volume supports the idea that they are related with low modes. And 
in fact it was shown in [23[ that the major difference between the best copy and the first copy lies in 
the low modes of the Faddcev-Popov operator which increase in number for a decreasing minimum of 
the gauge functional. 



3.7 Coulomb gauge results 

The propagators have also been measured on the lattice in the Coulomb gauge. The reason why it is 
specially interesting to consider this gauge is that relations between the behaviour of the propagators 
and the confining properties are particular ly tr ansparent. Let us first describe the continuum formalism, 
which was established by Christ and Lee jl3l| . The Yang-Mills theory can be written in terms of the 
Hamiltonian : 

H=\ J d 3 x (GUx)Gldx) + \ J d 3 xGt J (x)Gt J (x] 

+ [ d 3 xd 3 yp a (x)K ab (x,y)p b (y) (81) 



In this formula G± is the transverse part of the chromoelectric field (conjugate to A±), the kernel K 
can be expressed as the convolution of the Faddeev-Popov operator with itself and the colour density 
p can itself be written in terms of G± and Aj_ (in the quenched case which we are considering). It was 
further suggested by Gribov [36| that one should add to this Hamiltonian an extra term proportional 
to m 4 J d 3 xA±i(\7~ 2 )Aj_i in order to restrict the field configurations to what is nowadays known as the 
Gribov region. Doing so results in an equal-time transverse propagator D±(k, 0) oc |k|/(m 4 + (k 2 ) 2 .. 
At the same time, the time component -D44 of the propagator assumes the form 

D44 = V(x)S(t) + non instantaneous term (82) 

with V the vacuum expectation value of the k ernel K above. V(x) is different from the static quark 
potential, but, as has been shown by Zwanziger jl32| . it constitutes an upper bound for the latter. Thus 
a divergence of -D44 at large x is a necessary condition for the static quark potential to be confining. 

The Coulomb and Landau gauges can actually be considered as special cases of the so-called 
A— gauge specified by t he condition X djAj + (1 — A)9 4J 4 4 = 0. 

Some authors |l33l : Il34t 1 1351 : Il36l | have studied this general gauge and shown how the propagators 
pass continuously from the Landau to the Coulomb schemes. 

Since, however, those studies are still preliminary this section will be mainly dedicated to the case 
A = 1 (Coulomb gauge). 

The relationship between th e Lan dau and Coulomb gauge propagators has been studied in a more 



empirical way by Burgio et al [137J. Defining a Coulomb gauge scalar function £>c(|p|) as the po- 
intcgrated and rcnormalizcd scalar coefficient of the spatial tensor and comparing it to the Landau 
gauge function D(p 2 ) = G(p 2 )/p 2 they propose that the relation assumes the form 

\Pc(p)\ 

the UV asymptotic form of the rescaling function p — > pc (p) is known perturbatively. The important 
point, in what concerns the infrared region, is that pc{p)/p goes to a constant as p goes to zero. 
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Ref 


N 


year 




Lat. Size 


Di 


D 44 


ghost 


[140] 


2 


2000 


2.2 


14 4 ,16 4 ,.,,30 4 


.49-.51 


-1.9 


n.m. 


[63] 


2 


2003 


2.2, 2.3, 2.4,2.5 


16 4 - 24 4 


n.m. 


-2. 


n.m. 


[139] 


2 


2004 


2.2-2.8 


26 4 ,32 4 ,42 4 


0. 


-2.05 


-.245 


[141] 


2 


2007 


2.15-2.6 


36 4 


.12 


n.m. 


n.m. 


[138: 142; 143; 144] 


2 


2008 


2.15-2.6 


(24 4 ,32 4 ) 


.5 


-2. 


-1.22 


[145; 146] 


3 


2008 


5.9 


18 4 ,24 4 ,32 4 


n.m. 


-2 


n.m. 


Osaka [147] 


3 


2009 


5.8-6.2 


18 4 ,24 4 ,32 4 


n.m. 


-1.351 


-1.22. 


Berlin [1471 


3 


2009 


5.8-6.2 


18 4 ,24 4 ,32 4 


n.m. 


-1.13 


-1.22. 


[148; 149; 150] 


2 


2009 


2.2, 2.3, 2.4 


12 4 - 22 4 


n.m. 


-2. 


n.m. 


[151; 152; 153] 


3 


2009, 


5.8-6.2 


18 4 ,24 4 ,32 4 


.15 


-1.61 


n.m. 




3 


2011 


5.8-6.2 


Anisotropic 


.08 


-1.92 


n.m. 



Table 2 Summary of the infrared exponents of the gluon and ghost propagators from lattice simulations in 
Coulomb gauge; the notations are similar to the ones used in Landau gauge. N corresponds to the choice of 
gauge group. In refs. [1451 ; 11461] the results are given in x space rather than in p. The val ue "- 2" which we quote 
corresponds to the linear potential which they report. See in text for what regards rcfs. 139; 63; 148; 149; 150] 



This implies that a linearly vanishing Coulomb propagator, such as proposed by Gribov (see above), 
corresponds to an infrared finite and non- vanishing Landau gauge gluon propagator. 

Let us finally recall that the Coulomb gauge is not complete : even after the Coulomb condition 
diAi = has been imposed, it remains possible to perform a space-independent but time-dependent 
gauge transformation. Therefore a further step in gauge fixing is necessary; this residual gauge sym- 
metry is a continuous one ; taking care of it does not m ean that one does not have any longer to 
consider the discrete Gribov problem. It has been shown ( |l38j ) that this second step in gauge fixing 
is important in order to reduce scaling violations for the transverse propagator although it does not 
affect the ghost and temporal ones. 

3. 7. 1 Lattice results 

In table[2]we present the results of the several groups that have actually attempted to check numerically 
the expectations we have presented in the beginning of this section in a form similar to what we have 
done for the Landau gauge. The small p behaviour is parametrised as (p 2 ) a , meaning, in particular, 
that, for the transverse propagator, Gribov's suggestion corresponds to a = 0.5. Because of the very 
small number of simulations performed with SU (3) and since the general arguments about the IR 
behaviour appe ar to be independent of the gauge group we have included in this case the data obtained 
for SU(2). ref. |l39l ] measures directly the Coulomb potential rather than the temporal pr opag ator. 
Their result, an IR behaviour slightly more singular than p~ 4 may also have been seen in [138] but 
has not been reproduced in subsequent work on £>44. Greensite and his collaborators also measure 
the potentials (both coulombic and static interquark) in x-space. Their results agree with a linear 
long-distance behaviour and with the dominance of the Coulomb potential over the static one. 
A few words about the artefacts are also in order here. 

The scaling violations in the propagator, which are v ery i mportant as has been noticed for long 
and can be seen for example in figures 2 and 3 of ref. [153| . resist the usual cutoffing techniques 
devised to reduce the discretization effects. Their origin has been traced back to the definition of the 
instantaneous propagators, which induces a spurious dependence over |p|. A solution to overcome this 
difficulty consists in the use of anisotropic lattices in which the temporal spacing is much smaller 
than the spatial one. This is compensated for b y a la rger number of points in the time direction. The 
efficiency of his procedure has been checked in [153| . 

To summarize, the situation is not completely settled yet. However, in contrast with the situation 
regarding the Landau gauge, all Coulomb gauge simulations seem to be in qualitative agreement with 
Gribov's and Zwanzigcr's statements : 

• the transverse gluon propagator vanishes in the infrared. 

• the equal-time temporal propagator D44 diverges with R. 

Still, the precise infrared exponents are not exactly known yet. For D44 a linear divergence with 
R (1/p 4 in Fourier space) appears to be compatible with all data. Although there is a debate about 
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the value of the coefficient, it is admitted that it is larger than the known string tension. This is in 
agreement with the statement that -D44 is an upper bound for the static quark-antiquark potential. The 
crucial role of the l owest Faddev- Popov eigenvalues in building this confining potential is explicited very 
clearly in ref. [152| |. For the transverse propagator, t here is a problem with data giving an IR exponent 
smaller than 0.5. According to the argument of ref. }l37l ] this would correspond to a diverging Landau 
gauge gluon propagator, which is excluded. On the opposite no Coulomb gauge simulation results in a 
transverse gluon infrared exponent greater than .5, which, according to the same reasoning, means that 
all of them are in contradiction with any Landau gauge result predicting a vanishing gluon propagator. 
All these points still need to be clarified. 

On the other hand, as shall be discussed in sec. 14. f .31 the authors of ref. [3l| have very recently 
demonstrated that Gribov's formula for the equal-time spatial gluon propagator and the corresponding 
ghost dressing might be seen to admit bot h sca ling and decoupling behaviour for the resulting GPDSE 
in Coulomb gauge. Thus, as shown in ref. [l54j . the picture for both Coulomb and Landau gauge DSE 
solutions would be pretty the same, although the current lattice data, for the available momentum 
range, appear to be compatible with both classes of solutions. 

3.8 The dimension-two gluon condensate from the lattice 

As seen in sections 3.f and 3.2 and refs. [l55t Il56t If 57l : If 581 : [73|, the ghost and gluon propagators as 
well as the resulting strong coupling constant Eq. (fT3"|). do not run as perturbation theory requires 
in the energy range 2.5 - 7. GeV. This is surprising since it is widely believed that the perturbative 
regime is good above 2 GeV or at least 3 GeV. 

Having to deal with a non-perturbative correction to perturbative QCD we resort to the Operator 
Product Expansion (OPE) approach jl59j . In Landau gauge (which is the only one we shall consider 
in this section) there is only one dimension- two operator which has the vacuum quantum numbers: 
A% = A 2 . The fact that it is a dimension- two operator implies that it behaves as l/p 2 up to 
logarithms. This explains why the non-perturbative effects arc stronger when computed in Landau 
gauge (and more generally in any fixed gauge) than on gauge invariant quantities since the dominant 
gauge invariant gluonic operator is G^„ which has dimension- four and thus is cx I/(p 2 ) 2 up to 
logs. 

3.8.1 Several comments about OPE using A 2 

The use of OPE with A 2 in Landau gauge has been criticised. We would like to present our justification 
before going further. 

• A 2 is not a gauge invariant quantity although, in Landau gauge, it is invariant for infinitesimal 
gau ge tr ansformations as well as BRST transformations. As was already pointed by the authors of 
ref. 1 1 601 ] , it is legitimate to apply OPE with a gauge dependent quantity in a gauge theory. 
Indeed, to our knowledge, all the arguments used to prove OPE for a Lagrangian field theory can 
be applied to the theory defined by adding the gauge-fixing term to the gauge invariant QCD 
Lagrangian, including then the non-physical ghost fields to restore unitarity, i.e. to QCD in a fixed 
gauge. 

• {A 2 ) is ultraviolet divergent like the cut-off squared. As stressed in [l59j . when speaking of a 
condensate we think of the infrared modes of (A 2 ). How can wc discriminate in a theoretically 
sound manner the infrared modes from the ultraviolet ones ? A cut off in the loop momenta would 
be much too crude. The best is precisely to use the OPE expansion: given a quantity Q(p 2 ) we 
perform the following expansion 

Q(P 2 ) = Qpert(p 2 ,M 2 ) + C2 ilson (p 2 ,H 2 ){A 2 (n 2 )) + .... (83) 

A well defined renormalisation procedure and renormalisation scale is mandatory to be allowed to 
compare (A 2 ) computed from one quantity Q(p 2 ) and an other one Q'(p 2 ). Two different quan- 
tities will have different Wilson coefficients: Q pci t{p 2 , M 2 ) 7^ Q' pe rt(P 2 i M 2 ) an d Cwiison > ^ 2 ) 7^ 
C® ilson {p 2 , n 2 ), but the same (A 2 ). Indeed (A 2 ) is a property of the vacuum. It is not a lattice 
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artefacts, it is defined in the continuum limit, at vanishing lattice spacing. It depends on the gauge, 
and on the vacuum properties: the number and masses of the dynamical quarks. 

• The method just advocated consists in separating in a quantity Q(p 2 ) the perturbative contribution 
from the dominant non-perturbative one. However it is known that the perturbative series, i.e. the 
contribution of the identity operator in the OPE, is only asymptotically convergent. This means 
that higher order terms sum up in what is called "renormalons" which precisely behave like 1/p 2 
up to logs and are apparently non distinguishable form the effect of A 2 . It is not even clear that the 
distinction has a well defined theoretical meaning. This difficulty applies as well to gauge invariant 
operators an d th us to all the activity around what is called "QCD sum rules". This has been 
discussed in jl6l| . To make it short, the authors concentrate on the issue: can we use a condensate 
estimated from one quantity, say Q(p 2 ) for the expansion of another quantity, say Q'{p 2 ). They 
show that in Q(p 2 ) — cQ'(p 2 ) (c being a relevant prefactor) the renormalon ambiguity cancels. 
In other words, the properly defined difference between the two expansions contains a convergent 
perturbative series. If we stop at order n in the perturbative expansion, and assume that the sum 
of the expansion from n + 1 to oo is bounded by the term of order n, we get the following criterium 
of validity. We can indeed use the condensate estimated from Q(p 2 ) in the expansion of Q'(p 2 ) 
if the non-perturbative contribution (oc 1/p 2 ) is significantly larger than the last (highest order) 
perturbative contribution. This, of course depends on the energy. When the energy is very large, 
the last perturbative contribution dominates any 1/p 2 term and one can be satisfied with the 
perturbative series. In intermediate energies the l/p 2 -term dominates over the last perturbative 
one provided one has gone far enough in the perturbative expansion. At even lower energies all 
higher dimension operators contribute and OPE is no more applicable. 



3.8.2 Computing the Wilson coefficient 



A momentum dependent quantity Q(jp 2 ) with vacuum quantum numbers can be inverse power expanded 
as in Eq. (|83[) . Q can be the strong coupling constant, the quark field rcnormalisation constant, other 
renormalisation constants, etc |162| . The coefficients C® ilson are often called Wilson coefficients. All 

(p 2 , /i 2 ) are of the type 



Wilson coefficients CS-,„ 

wiison 
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(84) 
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From Eq. ((83)) and Eq. ((84)) we see that there is always the same factor (g 2 (fi 2 )A 2 (n 2 ))y^ s , (through- 
out this section we choose the MS scheme and the renormalisation scale fi = 10 GeV). We will therefore 
give the fitted values of the condensate as {g 2 {ii 2 )A 2 {)i 2 ))y^ s . 

In practice one can show that the best and most general fitting formula is: 
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At leading logarithm for the non-perturbative correction [15 
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for Nf = 2 
for N f =0. 



Notice that e is small. Therefore the corrective factor in Eq. (1871) is almost scale invariant. It is 
noticeable that the exponent is the same for all quantities Q(p 2 ) Ej- In Eq. (|87|l the only term which 
depends on the measured quantity is d^tree, which is given in table [3] for different quantities. 





gluon prop 


ghost prop 




£%3g sym 




3/32 


3/32 


9/32 


9/32 



Table 3 entree as defined in Eq. (|87|) . az g sym is the coupling constant extracted from the three gluon coupling 
with all momenta equal to p 2 [24| . 



Beyond the l eading logarithm, the Wilson coefficient have been computed at three loops by Chctyrkin 
and Maier jl63| for the propagators, and hence also for ar- 

3.8.3 Numerical results 

Once a quantity Q(p 2 ) has been computed on the lattice we need to treat the lattice artefatcs of order 
a 2 . There are two types of such artefacts, the ones which originate in the hypercubic geometry of the 
lattice, related to the group, and the ones which have the continuum geo metry, ( 4). 

The first type of artefacts are corrected via a non-perturbative method |ll7l : Ill9t Il64j j which fits 
from the data themselves for different orbits of the discrete group, the dependence of the considered 
quantity on the invariants of the discrete group which H4 are not invariants of 0(4). An extrapolation 
to the value zero of the latter invariants produces a result free of the hypercubic artefacts. Let us call 
the latter Q smoo th(p 2 )- 

Having got rid of hypercubic artefacts we still have to take into account the 0(4) invariant ones. 
The dominant one, at order a 2 is oc a 2 p 2 . We just add such a term to the fitting equation Eq. ([56"| : 



Qsmoothb 2 ) = Qpert(p 2 ,M 2 ) ( 1 + V)>5fc J + ^2 P 2 O 2 p 2 . (89) 

The function Q pm t{p 2 , M 2 ) is known from perturbation theory up to a multiplicative factor Qport(A* 2 , A* 2 )- 
The ratio C^ ilson /Q pe rt is also known, but (A 2 ) has to be fitted. Altogether, taking into account the 
coefficient c a 2 P 2, we have three quantities to be fitted from the function <5 sm ooth(p 2 ) which provides 
one data for each p 2 in our fitting interval. This is sufficient to perform a fit but there is a source of 
instability due to the fact that 1/p 2 decreases with p 2 , leading to a possible partial cancellation of the 
contribution of the artefact and the condensate. However the lattice spacing dependence of the term 
oc a 2 p 2 and the condensate cx 1/p 2 up to logs are very different, which allows to check the quality of 
the fit. 

We will now consider several Q(p 2 ) and provide the estimated values of the condensate. This is 
shown in Tab. 0) For the particular case of cxt, which we paid special interest to in this work, we also 
borrow from ref. (73| the plots shown in Fig. [5J where the very good agreement between the prediction 
from Eqs. (|86I87|) and the lattice data for the Taylor running coupling is manifest. 

As will be also mentioned in App.|21 this dimension-two gluon condensate is seen to have an impact 
on the ghost-gluon vertex and has to be taken i nto a ccount when solving the GPDSE and comparing 
with lattice data for the ghost dressing function [l65j . 

15 For the coup ling constant co mpu ted from the three gluon coupling with one vanishing momentum, the 
formula Eq. (JHTJ) does not apply |157lj . 
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fitted g 2 R (A') 



order g 2 {A 2 ) 


gluon propagator 


ghost propagator 


CtT 


3- gluon a asym 


3-gluon a sym 


Tree 


2.7(4) 


2.7(2) 








LL 






5.2(1.1) 


10(3) 


6.8(1.5) 


0(a 4 ) 






3.7(8) 







Table 4 Comparison of estimates of g 2 (A 2 } from different quantities. All are taken at th e scale /i — 10 GeV. 
Tree means tree level for the Wilson coefficient. The data in this line come from [158]. LL means leading 
logarithm for the Wilson coefficient. The data in this line are taken from ref. [73^ . 0(a 4 ) refers to Chetyrkin 
and Maier computation. 
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Fig. 5 (Left) Comparison, borrowed from ref. [TJi °f lattice estimates for the running coupling in Taylor 
scheme and the non-perturbative prediction including the dimension- two gluon condensate, (A 2 ); qt in vertical 
axis is plotted in terms of the square of the momentum in horizontal axis. (Right) The same at very high 
momenta to reach the perturbative regime, where the gluon condensate effects can properly be neglected, and 
check the consistency of the results, in particular the estimate of A^g. 



4 DSE numerical solutions 

4.1 Solving the ghost propagator DSE 

Let us analyze first in this section the general picture for the low-momentum Landau-gauge DSE 
solutions by applying the approach proposed for the first time in ref. [2(1]: to combine lattice gluon 
propagator results with a ghost propagator DSE truncated by a well supported hypothesis for the ghost- 
gluon vertex. The same procedure has been very recently followed for the analysis of the Coulomb-gauge 
DSE solutions in ref. |31| leading as w e shall discuss below, to a very analogous picture. It should be 
noted too that the authors of ref. |l66l | also successfully applied the same strategy of combining lattice 
data and DSE, in particular by invoking the lattice data for the truncation of the gap equation and 
studying the chiral symmetry breaking. 



4-1.1 Landau gauge solutions 

We shall present now the results of ref. [26| for the ghost solutions in Landau gauge. The goal was 
to see whether the two types of solutions (of = and lap + etc = 0) suggested by the previous 
analytical discussion in section [2] actually exist for the same gluon propagator input. A positive answer 
came out by solving numerically the ghost SD equation for given gluon propagator and vertex, as we 
shall briefly describe in the following. To this goal, we invoke again the renormalized DSE given by 
(HHJ) and cast it into the appropriate subtracted form: 
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F R (k 2 ) 



G R ((q 



(2tt) 4 
-kf) 



£*((<? -fc') 2 ) 



F R (q 2 



(90) 



where we work in the MOM scheme, and set k 1 appearing in ea. (jl9p as the squared renormalisation 
scale /i 2 (fi has been chosen at an optimum 1.5 GeV, not too high to allow the lattice data to be safe, 
and not too small in order that the differences between solutions at small momenta can be clearly 
displayed). An IR finite gluon propagator (olg = 1) extracted from lattice data in pure Yang-Mills 
theory, with Wilson gauge action, j3 = 5.8 and a lattice volume equal to 32 4 , is used for momenta lower 
than 1.5 GeV ; this choice is justified to have moderate UV artefacts. This is then extended to larger 
momenta using a one loop asymptotic expansion (with Amom = 1 GeV corresponding to the standard 
Ajj s = 0.240 GeV of lattice quenched QCD). As for the ghost-gluon transverse form factor in Eq. (flSl) . 
Hi(q, k), it is taken to be constant with respect to both momentaF^I. As we said above, this is suggested 
by the lattice data for q = k (i.e. for zero gluon momentum), but we extend it to all values of q and k. 
The authors of ref. |109| find a bare vertex very close to 1 in this zero momentum gluon configuration 
for a large range of \/cp '. Note that we have re-defined the coupling as g 2 = NcZig^fi 2 ) in Eq. (1901) . 
which, as far as the constancy for the ghost-gluon transverse form factor is assumed, depends only on 
the renormalisation point chosen for the propagators; it is furthermore independent of the particular 
way used to define the renormalisation of the vertex. On the other hand, the redefined coupling g can 
be also written in terms of bare quantities or related to the well-known running coupling in the Taylor 
scheme: 



g 2 = N c g 2 R Z x H 1R = N c g 2 B Z 3 Z 2 /Z 1 H 1R 
= N c g 2 B ZzZl Hib 

= N c &(ji 2 ) Hib (91) 

where it should be remembered that H\ R = 1, in MOM scheme. Eq. (jHOJ) can be still transformed 
to a new form which makes the numerical calculation and the presentation of the various solutions 
easier; for this, we subtract the equation at k = 0, to let the value of F R (k) at the origin appear and 
to eliminate the reference to the particular renormalisation point fx, and we redefine also the unknown 
function to be calculated as F(k) = gF R (k). Then the reference to the value of g also disappears; we 
end with : 



d 4 q 

F{k 2 ) F(0) J (2tt) 4 
G R ((q-k) 2 ) 
((q~k) 2 ) 2 



1 jk-q)'- 
k 2 q 2 

G R {{q?) 
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F(q 2 ) 



(92) 



Equation Eq. can be solved for F(k 2 ), for a set of values of F(0). It is easy to see that, from 
this, the desired solution of eq. (|90| can be straightforwardly reconstructed for any renormalisation 
point and any value of g. Indeed, as the MOM rcnormalization condition imposes g(fi) — F(fi 2 ), for 
any given [i and <7, we have just to identify the value of F(0) such that F(fj, 2 ) = g and reconstruct 
then F R (k 2 ) through 



F R {k 2 ) 



F(k 2 ) 

<?0) 



(93) 



16 This cannot be an exact statement, as already shown in perturbation by the calculations of ref. [4|; [78|] 
although finite, the vertex invariants do depend on the momenta through the running a s . 
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By construction, all the solutions found in this way are finite at the origin and they correspond to the 
"decoupling" family of solutions that were described above in section[21 On the other hand, the solution 
which diverges at vanishing momentum appears as an end-point for the solutions of this "decoupling" 
family that can be approached by making F(Q) larger and larger so as to get the limiting case: 0. 

F(0) 

Then, the critical "scaling" solution will be found by setting -J— = in cq. (j9"2")l . As we will discuss 

below, this strategy for solving the DSE, after it had been followed in [26j], was further applied by the 
authors of ref. [28| to the analysis of the coupled ghost and gluon DSE in Landau gauge, where they 
also concluded that the choice of F(0) amounted to fix a boundary condition for the DSE system and, 
consequently, to determine the class to which the actual solution belongs : either a "decoupling" one 
for any finite value of F(0) or the unique "scaling" one for F(0) — > oo. Unfortunately, the authors 
of ref. [28| missed the connection of F(0) and g, and hence with the coupling at the renormalization 
point, gn(u?\ given by Eq. ([93]). 

In ref. [261 ]. the solutions of ea. ([9"2")) with the integral cut in the UV at q = 30 GeV were obtained 
after discretization in k and q and solving by iteration^. The results are the following: 

1) Critical case, scaling solution: One finds a solution with -^i— = 0, i.e. F(Q) = oo, the 
corresponding "critical" constant being 

g c 2 = F(1.5 GeV) = 33.198.... (94) 

The relation of ea. (|43|) for ag = 1, obtained in sec. [2 should be verified by the numerical solutions 
and it happens to be very well satisfied: 

£ 2 ,4V )G (2) (0) 1 ^ 0.994.... (95) 

The integration near k = can be improved by taking explicitly into account the analytical behavior 
of the kernel, and assuming that the solution behaves as 1/k at small k. This imposes eq. (|4"3")) , and 
one indeed can check that 

l im n T7^^ 2 ^ 2 ) 2 G (2) (fc 2 ) 06) 

although very slowly. 

2) Regular case, decoupling solution: One finds a solution and only one for any F(0) > 
with the method of solution described above. As can be seen below, a numerical solution at g 2 ~ 29 
corresponds to the best description of lattice data (see Fig. [7J borrowed from ref. HH). Furthermore, 
the asymptotic low-momentum behaviour for t he g host dressing function, giving a next-to-leading 
£; 2 log(£; 2 ) term, (cf. Eq. (j5Tj) ) is checked in ref. [26[ (this is shown in FigHJl also borrowed from this 
work) and the slope appears to agree pretty well with the prediction from Eq. (f5"Tj) : 4.06 against 4.11. 



The critical value of the coupling constant, as well as the corresponding curve of F(k), can be very 
well approximated by the regular decoupling solutions at very large F (0). When F(0) is larger and 
larger, Eq. (|5ip remains valid only in a smaller and smaller region near k = while, in an intermediate 
region, one observes the expected 1/k behaviour. In this way it is possible to show that, as long as the 
coupling constant does not exceed a critical value g C rit, where the scaling solution emerges (see also 
sec. 14.1.41 ). the solution goes to a constant in the infrared : one is in the decoupling case. As soon as 
the coupling constant reaches that critical value the solution converts to the infrared-infinite one 
("scaling situation"). Actually, the DSE solution accounted for the scaling solution when g cr it — 33.2 
at the renormalization point fi = 1.5 GeV, as it is shown in Fig. [7] Admittedly, a real resolution of the 
DSE equations would require solving reciprocally the DSE for the gluon and verifying the compatibility 

17 It should be noted that minus the integral in the r.h.s. is positive, allowing an easy convergence. We linearize 
it at each step, following the Newton method, to accelerate the convergence of the iteration procedure, as 
suggested by Bloch [68| . 

18 The existence and the value of this critical coupling had already been noticed, see for instance ref. |68lj . 
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Fig. 6 The a + bk 2 log(fc 2 ) fit at small momentum (dashed line) to our continuum SD prediction for the ghost 
dressing function, renormalised at fi — 1.5 GeV for g 2 = 29. (solid line) ; the slope of the fc 2 log(fc 2 ) term is 
4.06 ; the agreement with the expected coefficient of fc 2 log(fc 2 ), 4.11 from the eq. (|51[) . is striking. 



of the solution with the gluon dressing function which was used as an input to build the kernel for the 
same g- value. We will deal with this below. 

In conclusion, in the case ag = 1, a continuum set of IR finite decoupling solutions for arbitrary 
F(0) emerges, and a unique singular scaling solution for g 2 = g 2 , with ap = — A, which appears to be 
the end-point for the previous ones. 

4-1-2 Comparison with ghost lattice data 

After the analytic study of the general solutions of the GPDSE in sec. [2] and the previous numerical 
analysis exhibiting both types of solutions for the ghost dressing functions, either regular (decoupling) 
or singular (scaling), the question also addressed in ref. [26j is: which one is effectively realised on the 
lattice, and therefore in true QCD? 

A better means to provide us with an answer is offered by the numerical calculation in previous 
subsection that, predicting the behavior of the respective solutions for the ghost over the whole range 
of momentum, can be confronted to the lattice estimate for the ghost dressing function and that can 
lead us to identify which one offers the best agreement with the data0. This is done in ref. [26| and 
displayed here in fig. where one can see that the singular scaling solution appears to be clearly 
discarded by the lattice data around and below k = 0.5 GeV. On the contrary, a very good description 
of the lattice data in the range g 2 = 28.3 — 29.8 (clearly below the critical value) is found. This striking 
agreement, although for a narrow momentum window, is illustrated by Fig. [7] (As an indication, we 
quote the IR limit Fr(0) = 2.51 for the same fx = 1.5 and g 2 = 29). 

An additional consistency test is obtained from using Eq. (|91[) to connect the continuum g 2 to the 
lattice bare quantities, 

g 2 = N c g 2 R z x = N c ~F 2 ((x 2 ) G B (/i 2 ) H 1B , (97) 



19 At this stage, it is useful to stress the advantage of working with the renormalised form of the SD equations; 
inde ed the continuum and lattice versions are more directly comparable than the bare ones. As the authors of 
ref. [24| have shown, the bare lattice equation for the ghost is affected by an important artefact which vanishes 
only very slowly with the cutoff, being of order 0(g 2 ). In the renormalised version, this effect is included in 
the renormalisation constant Z3, and we are left only with the much smaller cutoff effects of the type 0(a n ). 
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dressing function lattice result 
best SD solution 
singular SD solution 




Fig. 7 Comparison, borrowed from ref. [2(||, between the lattice SU(3) data at j3 = 5.8 and with a volume 32 4 
for the ghost dressing function and our continuum SD prediction renormalised at /i = 1.5 GeV for g 2 = 29. 
(solid line) ; the agreement is striking; also shown is the singular solution which exists only at g 2 = 33.198.... 
(broken line), and which is obviously excluded. 



and then checking whether our range g 2 = 28.3 — 29.8 is reasonably consistent with the r.h.s. of ea. (l971) 
as given by lattice data. In spite of the crude approximation made to be left with Eq. (f9"T|) , the 
result of the check is very encouraging: indeed, from the above value of g 2 found in the continuum on 
the one hand and the lattice data (3 = 5.8, Gb{^ 2 ) — 2.89 and Fb{p?) — 1-64 is here chosen as 
1.5 GeV) on the other, one finds H\b — 1.2 to satisfy equation (|9~71) . which represents some kind of 
average on momenta. This number should be compared to the SU(3) lattice measurements which are 
performed at q = k and have been submitt ed to a renormalization such that the result takes the value 
1 at q = 3 GeV, and which give about 1.1 |l09j |. with large errors. 

Another striking way of presenting the difference between the regular decoupling solution and the 
singular scaling one is in terms of the familiar product, G i?,(fc) F^{k) 2 , which is proportional, at least 
for the UV dom ain, to the running coupling in Taylor scheme. This is shown in Fig. |SJ 

In ref. 105j new ghost propagator data are provided, coming from larger volume lattice simula- 
tions and covering a wider momentum range, with smaller momentum data, to compare with. The 
confrontation of those data with the results from the integration of the GPDSE truncated with the 
help of a lattice-based gluon propagator as done in sec. 14.1.11 requires, to account succesfully for the 
very low-momentum data, to go beyond the hypo thesis of constancy of the ghost-gluon transverse 
form factor, H\. A possibility, studied in ref. [165| consists in applying an ansatz for the transverse 
form factor, H\, which is inspired by the OPE description for the non-perturbative corrections of the 
ghost-gluon vertex, obtained by applying the same procedure outlined in sec. 13.81 and by taking the 
gluon condensate value of table [4] 



4-. 1.3 Coulomb gauge solutions 



As we previously 
the one in ref. [2J 



mentioned, the authors of ref. [311 ] recently performed a study (very analogous to 
| for Landau gauge) of the GPDSE in Coulomb gauge, obtained within the (second 

20 First, it is valid up to finite cutoff effects, as well as volume effects; second, we have replaced the lattice 
vertex invariant HiB(q,k) by the constant H\b which, as we discussed above, is a rough approximation over 
the momenta which are actually implied in our calculation. In app. [X] some non-perturbative c orrec tions for 
this invariant, that appear to be pretty well in agreement with some SU(2) lattice estimates |165lj . will be 
discussed. 
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k 2 [GeV 2 ] 



Fig. 8 The same than in Fig. [7] but for Git(k)Fii(k) 2 ; up to a factor g 2 /4ir, the scaling (dotted) curve 
corresponds to the a(k 2 ) presented in Fig. 8 by Fischer [1671] : the shape is very similar. 



order) functional formalism, in order to investigate the low-momentum ghost dressing solutions. They 
took Gribov's equal-time spatial gluon propagator dressing function PH. 



dk 4 G(kl,k 2 



&4 



2 7k 3 



(98) 



as the input required to build a kernel and solve the GPDSE, again with the approximation of replacing 
the fully dressed spatial ghost-gluon vertex by the bare one (this is, also in Coul omb gauge, an exact 
result in the limit of a vanishing incoming ghost up to all perturbative orders [l69l |). Thus, the GPDSE 
is rewritten as follows: 



9M 

F(k 2 ,fi 2 ) F(p 2 ,(i 2 ) " u (4tt) 2 
where / represents the angular integration, 



r^F(q 2 , M 2 ) (/(k 2 ,q 2 ;m)-/(p 2 ,q 2 ;m)) , (99) 
■Jo Q 



/(k 2 ,q 2 ;m) 



dz I 




-1/2 




-1/2 



(100) 



It should be emphasized that the ghost propagator dressing function in Coulomb gauge is st rictly 
independent of the energy, fcf , as a non-perturbative consequence of the Slavnov- Taylor identities [l70j . 

Assuming a pure powerlaw behaviour, -F(k 2 ) ~ (k 2 ) QF , for the ghost dressing function and ana- 
lyzing asymptotically Eq. (|99)l. one is left in Coulomb gauge again with the two same cases we have 
encountered in Landau gauge: (i) ctF — (decoupling) and (ii) ap (scaling). As well in Landau 
as in Coulomb gauge, a massive gluon propagator generated via the Schwinger mechanism or Gribov's 
fomula for the equal-time spatial dressing leads to ac = 1 and thus olf = —1/2- In particular, from 
the next-to-leading analysis in sec. 12.21 of Eq. ((24|) . one obtains [l?} : 

21 In ve ry good agreement with the Euclidean SU(2) lattice results obtained for small lattice couplings in 
ref. [T6l. 



37 



IOtt 2 



NcH igR (ji*)B0j* 



1/2 



1/2 



1 



aT(0) M 2 



In ■ 



A'/ 2 



11 

IT 



(9 



M 4 



if a F ^ , 



if = 



(101) 



If cvi? ^ 0, the perturbative strong coupling defined in the Taylor scheme [72 
a non-zero constant at zero-momentum, 



a T = fl , |-/(47r), goes to 



lim axiq 2 

q 2 ->-0 



lim 



9 



2 (M 2 ) 



47T 



2 z% 2 V)^W) 



57T 



2iY c #i 



(102) 



as can be obtained from Eqs. (|54I101[) . In the case = 0, the subleading correction to the non- 
zero finite value for the zero-momentum ghost dressing function, given by Eq. (|101[) . is controlled 
by the well-defined zero-momentum limit of c¥t(<Z 2 ) = (M 2 /q 2 )c t r(q 2 ), which is the extension to the 
Taylor ghost-gluon coupling case [2!| of the non-perturbative effective charge defined from the gluon 
propagator in ref. (74| . 

The same two cases result from the analysis of Eq. for the Coulomb gauge in ref. [3l[ , where a 
ghost propagator dressing function behaving asymptotically as either a constant or F(k ) ~ (k 2 ) -1 / 2 
is analytically found and confirmed by a numerical study. This can be seen in the left plot of Fig. [5] 
which we borrowed from ref. [3l| . In Coulomb gauge, the lattice results for the ghost propagator, within 
the available momentum window, may agree with decoupling and scaling solutions, and cannot help 
to discriminate. 

It should be noted that the lagrangian approach to the Coulomb gauge (the continuum functional 
formalism is based on the QCD Lagrange density) is not the most widely used. However, and despite 
the technical difficulties mostly related to the inherent non-covariance o f the Coulomb gauge, some 
recent progresses have been made in order to derive explicitely the DSE 169t_ 170l ; Il7lj (allowing the 
previous analysis) or studying the Bethe-Salpeter equation for heavy quarks [172l ] . The method which 
is most widely applie d in the continuum is the canonical formalism based on the QCD Hamiltonian 
density operator [32|; Il73t Il74t 1 1 T5t Il76l |. Dyson-Schwinger-like equations (in a space with one less 
dimension) for the equal-time correlators are obtained in the canonical formalism and again the two 
types of solutions, critical and subcritical, are found for the propagator dressing function [32| . It 
seems to happen that, in the Coulomb- gauge canonical formalism, two different values for the infrared 
exponents emerge in the critical case and that the favoured one is the one which produces the most 
singular ghost dressing, which diverges as l/|k| in the infrared, i.e. similarly to the critical solution 
discussed above (with = —1/2). 



4-. 1.4 The critical limit of decoupling solutions from the GPDSE analysis 

In summary, the GPDSE in Eq. (f!M)) with the input of a gluon propagator borrowed from lattice QCD 
calculations can be numerically solved and two kinds of solutions result and appear to be controlled 
by the size of the coupling at the renormalization point Q git 1 )- For any coupling below some critical 
value, g C rit, an infinite number of regular or decoupling solutions for the ghost dressing, behaving 
as Eq. poip indicates, are found; for <?(/i) = <7 C rit, a unique critical or scaling solution behaving as 
Eq. (jlOip is found, and no other type of solutions appears to exist. In ref. [26j | . for a subtraction point 
/i = 1.5 GeV, a critical coupling g cr ;t ~ 3.33 and a very good description of ghost propagator lattice 
data with a regular solution of Eq. ((24)) for g(/j,) ~ 3.11 were obtained. 

As we shall discuss in the next section, these results were also recently confir med I l77l| by studying 
the coupled system of ghost and gluon propagator DSE in the PT-BFM scheme |l78j . This last work 
paid attention to the critical solution limit by studying how _F(0,^t 2 ) d iverg es as g(fJ>) — t <?crit — 

1.51, 

with a subtraction point /j, = 10 GeV. In addition, the author of ref. [154| applied the perturbative 
definition of the Taylor strong coupling in Eq. (|102D to compute this coupling with the gluon and ghost 



In QCD, one needs to provide a physical scale and a standard manner to proceed is by fixing the size of 
the coupling at a given momentum scale. This can be seen as a boundary condition to solve the DSEs. 
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Fig. 9 (Left) Ghost dressing in Coulomb gauge plotted for different values of the inverse of boundary condition 
fixed for the dressing at zero- momentum, -T(O) in terms of x = k 2 /m 2 , where m is the Gribov mass in Eq. <198|) , 
This plot is borrowed from ref. [3j . (Right) The Taylor coupling defined by Eq. (|102[) and computed from the 
gluon and ghost propagator results of ref. [H3 for g 2 : (/x) / '(4tt) = 0.16, 0.17, 0.18, 0.1812, 0.1815, 0.1817, with the 
subtraction point fi = 10 GeV. The curve for the critical limit is obtained by applying the results of ref. [261] . 
as explained in the text. 



solutions of 177f in order to s ee h ow the critical limit is approached. This is shown by the right plot 



of Fig. [51 borrowed from ref. 154| , where it can be also seen that all the curves for ax obtained for 
different values of <?(/i) tend to superpose over each other as g 2 //x 2 increases (right). As a striking check 
of consistency, the curve for the critical limit in this right plot of Fig. [5] is obtained by rescaling, up 
to giving Q!t(0) from Eq. (|102p with H± — 1 at zero- momentum, the results at the critical limit for 
q 2 A(q 2 )F 2 (q 2 ) numerically obtained in ref. [26| and plotted here in Fig. |5]of the previous subsection. 
Indeed, the critical value for the coupling at [i = 10 GeV can be read fro m th e critical curve in Fig. [9] 
and one gets g{fj) ~ 1.56, in fairly good agreement with the value of ref. [l77j . 

On the other hand, the numerical analysis of Eq. (|99| for the Coulomb gauge in ref. [3l[ also 
shows both the regular and the critical solutions to exist, but being controlled by F(0, n) (or r(0, /i) = 
1/.F(0, /i)) as a boundary condition with the size of the coupling fixed to be g 2 (n) = ~g 2 = An x 0.1187 
for Nc = 3 (see the right plot of Fig. [5]). Regular or decoupling solutions appear for finite values 
of F(0,fi 2 ) and critical or scaling solution for F(0,/i 2 ) — > oo. Nevertheless, the authors of ref. [|l| 
concluded that, as far as all the solutions join each other in the perturbative domain (see the right 
plot of Fig. IH]) and can be found for a fixed coupling, the boundary condition is not connected to the 
renormalization and claimed for a contradiction with the Landau-gauge results of ref. [26j . 

Finally, the two pictures for the low- momentum solutions from GPDSE, either a family of Landau- 
gauge DSE solutions labelled by the size of the coupling at the renormalization point or a family of 
Coulomb-gauge ones labelled by the zero-momentum ghost dressing val ue as a boundary condition 
independent of the renormalization, were reconciled by the work of ref. [l54j . The key point stems 
from the different renormalization prescriptions applied to t he gh ost propagator in both analyses. 
MOM scheme in the Landau-gauge analysis of refs. [2(| H?! Il77l |. and the prescription applied to 
the ghost propagator in eq. (3.20) of ref. [31| for the renormalization constant Z C (A, \g,r(0)]), where 
r(0) = 1/F(0, (x 2 ). In particular, this last renormalization constant depends on the boundary condition, 
-T(O), in such a manner that the value for this boundary condition is rescaling the ghost dressing 
function (and, as can be clearly seen in Fig. 2 of [3l|, it does not take the t ree-le vel value, 1, as 
happens in MOM prescription for the subtraction point). Thus, according to ref. [154| . the non-trivial 
connection between solutions in both schemes, that relies on the relation previously shown by Eq. (|93l) , 
comes out from the following property of Eqs. (|24I99|) : let F(q 2 7 fi 2 ) be a MOM solution of Eq. ((24)) 
for arbitrary coupling, g(fi); if we then apply the following transformation: 



g(jjL) -> s g(fi) , F(q 2 , v 2 ) -> -F(q 2 7f i 2 ) 



(103) 
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for any c-number s, the transformed dressing function verifies the DSE equation with the transformed 
coupling (of course, MOM prescription implies s = 1). Then, one only needs to choose s = g~/g(fi) and 
to apply the transformation to every solution of the MOM family and one will be left with a one-to-one 
correspondence between these solutions and the new ones 

F(q 2 ,^ 2 ) ee g M F (q 2 ^ 2 ), (104) 
9 

for the fixed coupling g, which can be identified by the zero-momentum value, -F(0,^ 2 ). This new 
family of transformed solutions obeys the same pattern as the Coulomb gauge family in ref. [3lj and 
corresponds, up to the fixed number g (which does not even depend on the renormalization point), to 
the family of solutions F(q 2 ) obtained by the analysis for Landau gauge in previous sec. 14.1.11 It is 
interesting to note that the strong coupling defined in the Taylor scheme can be also obtained from 
the transformed solutions as 

aria 2 ) ee t- q 2 D(q 2 , fi 2 )F(q 2 , fi 2 ) ee 9 -M q 2 D(q 2 ^ 2 )F(q 2 , M 2 ) , (105) 

47T 47T 

although it is obvious that neither F nor the coupling are in MOM scheme. 

Thus, the same picture for the low- momentum Green function solutions emerges in both Landau 
and Coulomb gauge from the analysis of the GPDSE: a family of MOM-rcnormalizcd regular decoupling 
solutions, characterised by the value of the coupling at the renormalization point; and a singular scaling 
solution as an end-point for the fa mily of regular ones. An interesting final remark is that the input 
parameter for the solutions in ref. [3l|, the zero- momentum ghost dressing, can be put in connection 
with the Gribov problem [3l[; while, for Landau gauge and MOM scheme, g(fi) is related to the strong 
coupling in Taylor scheme, as was shown in Eq. (|9T|) . This is not the case for the size of the fixed 
coupling, g, after applying Eq. (|103[) which is physically meaningless. 



4.2 The ghost and gluon propagator coupled DSEs 

Satisfying the GPDSE, as was required in sec. [3J is a necessary but not sufficient condition for a DSE 
solution to exist. Of course, the existence of a solution can only be confirmed by treating the infinite 
tower of DSEs, but this is an impossible task. In the previous section, sec. 14. 1[ this infinite tower of 
DSEs was truncated by plugging into the one among them to be solved, the GPDSE in that case, 
the available lattice data, or a model compatible with them, for the gluon propagator and the ghost- 
gluon vertex. On the other hand, the usual approach consists in applying a truncation scheme based 
on hypotheses and approximations that preserve the main properties of the theory and that leave us 
with a closed system of equations to deal with. The former approach can be seen to provide with a 
consistency analysis of the lattice and DSE picture for the solutions and benefits of not "polluting" 
the conclusions with the possible implications of any particular truncation scheme. However, the DSE 
picture should be completed by also applying the latter usual approach. 

As we shall discuss below, both scaling and decoupling solutions have also been proven to emerge 
when the DSEs are truncated so as to give a coupled system for the ghost and gluon propagators. 

4-2.1 Scaling solutions 

As we have repeated insistantly in this paper, it has been recognized for a long time that the set 
of solutions of the DSEs for the ghost and gluon propagators consists in a continuum of so-called 
"decoupling" solutions augmented with a unique "scaling" one (cf. ref. [66|); which one is actually 
encountered depends on the value of the coupling constant. Nevertheless, for quite a time, attention has 
mainly be paid to the scaling one which, in practice, was obtained by replacing the fully dressed vertices 
by ansatze which take into account as much information as possible (see ref. [ll| for a first review on 
the subject). The lo ops in DS Es had been proven to be dominated by the infrared contributions for 
the scaling solution jl08l : Il79j which had been fully worked-out in ref. [66[ (see also ref. [1^]). The 
infrared exponents for the power behaviour on the momentum for both ghost and gluon dressing 
functions being related by 2ctF + <^g = 0, the value for one of them, usually the one for the ghost, aj?, 
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completely characterizes the low-momentum solution. Under the assumption of a constant ghost-gluon 
transverse form factor, the only solution to emerge in t he interval [—1, —1/2] ( [66]) is ap ~ 0.595, as 
was numerically put forward by the authors of ref . |l80j , and independently confirmed by the analysis 
performed in ref. [5(| with the help of renormalization group methods (RGE). Then, the uniqueness 
of the the above-mentioned low-momentum solution defined by ~ 0.595 was discussed in two 
papers [10; [EH > nrst (wrongly) claimed to be true in general for the coupled DSE system [3(| and later 
on concluded to happen only for scaling- type solutions [HJ, i-e. provided that the relation 2ap + ac = 
is "a priori' assumed (we will pay attention to this in a next subsection). 

Very recently, the authors of ref. [28| re-analysed the problem of the low-momentum properties 
of the Yang-Mills Green functions by following both DSEs and RGE approaches and found both 
scaling and decoupling solutions to exist. They paid special attention to the truncation schemes and 
also claimed that only the (unique) scaling solution satisfies BRST invariance, while the decoupling 
ones would be at odds with it. This was presented as an incitation to prefer the scaling solution as 
the "real" QCD one but, as was previously mentioned when discussing the Gribov-copies problem, 
Gribov or Gribov-Zwanziger (either refined or not) approaches to avoid the copies already imply a 
BRST breaking and this only prevents the Kugo-Ojima confinement scenario from working. Other 
confinement scenarios are of course possible and nothing indeed prevents the decoupling solution from 
being, as lattice appears to indicate, the real QCD one. 

The properties and implications of the scaling-type low-momentum solutions have been extens ively 
discussed in the literature. We address the interested reader to reviews like the ones in ref. [TBI; EH; Il67| 



as well as to the numerous works qoted a bove or to others like ref. [1811 ] focusing on the Kugo-Ojima 
criterium implications, refs. Il75 :| l82tll83j ] on the infrared behaviour of vertices, ref. [8lJ on the analysis 
of infrared sigularities or ref. [184j about the study of scaling solutions in the maximally abelian gauge. 

We will now end this section by adding a few words about the elusiveness of decoupling solutions, 
after the scaling one has been proposed. 



4-2.2 Why have the decoupling solutions been so elusive? 



The decoupling or regular solutions have been missed for almost ten years. Why? they could not have 
been previously obtained by the proponents of the relation (Rela) because, as it seems to us, they 
discarded them from the very beginning, and thereby chose the critical value of the coupling constant, 
by making an implicit assumption when solvin g the so-called "infrared equation" for the ghost SD 
equation. This can be seen, for instance, in ref. [l79j . eqs. (43) and (44), or in the detailed discussion 
of Bloch eqs. (55) to (58). 

Let us explain this briefly. They consider the above unsubtractcd equation (note that this requires 
then an UV cutoff, which we avoid in our previous analysis by considering the subtracted form); we 
write again the unsubtractcd form: 



1 



F R (k 2 



= Z 3 - N c g R zi 



d 4 q 



1 - 



G R ((q-k) 2 )H 1R (q,k) 



k 2 q 2 



F R {q 2 ) 



(106) 



One must try to match the small k 2 behaviour of the two sides of Eq. (|106[) . This is done for example 
in eq. (58) of [68[- A condition is then written which consists in equating the coefficient of (k 2 )~ ap 
with the corresponding one in the r.h.s.. However, one notices that on the r.h.s., there is a constant 
contribution oc (fc 2 ) . Therefore unless the constant term Z% is cancelled by the integral contribution 
for k — > 0, we have necessarily ap = 0. To have ap < as the author finds, one needs this cancellation. 
This is what is implicitly assumed, but not stated explicitly. The condition of cancellation is : 



Z3 = Ncg^zx 



d 4 q 
(2tt)4 



1 - 



(k.qf 
k 2 q 2 



F R (q 2 ) 



(107) 



k=0 



However, this additional equation does not derive from the starting SD ghost equation, and indeed it 
is not satisfied in general by the solutions of this basic equation, as we show by displaying actually IR 
finite solutions. In fact, it can be valid only for a particular value of the coupling constant, the 
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critical one which is solution to the equation of Bloch, his eq. (58), and which we derive rigorously 
through the subtracted equation. A similar conclusion is obtained in the analysis of ref. [28[ , although 
its authors missed the connection between their boundary condition, the zero-momentum value of the 
renormalized ghost dressing function, and the coupling size at the renormalization momentum. 

4-2.3 Decoupling or massive solutions 



A decoupling behaviour has also been proven to result as a solution of a coupled system of gluon 
and ghost propagators DSEs. First, the authors of ref. |185| implemented some ansatzc based on 
Slavnov- Taylor identities for the involved full vertices, applied a particular angular approximation 
when integrating the ghost self-energy and thus obtained a "massive" gluon propagator {olq = 1), 
although they claimed this to be compatible with an enhanced ghost propagator. Then, as mentioned 
above, the Schwinger mechanism of mass generation [6^] was proven to be consistently incorporated 
into the gluon propagator DSE through the fully-dressed non-perturbative three-gluon vertex and to 
give rise to the generation of a dynamical gluon mass ffdt li~9l. Then, a massive solution for the coupled 
ghos t and gluon propagator DSE was shown to apppear [71j in the PT-BFM truncation scheme (see 
also |l86j ). The lattice data result to be furthermore very well accommodated within coupled DSEs 
in the PT-BFM scheme [20l|. As a matter of fact, as will be seen below, the PT-BFM DSEs solutions 
have been shown to asymptotically behave as Eqs. ([54l [55)1 predict for a decoupling solution 177 1. 
The authors of ref. [28| also confirmed the decoupling solutions to be present by the analysis of the 
coupled DSEs and the functional Renormalization group equations (FRGs). They also obtained an 
infinite family of decoupling solutions, as was discussed in sec. 14.1 .41 but they used the zero-momentum 
ghost propagator as the boundary condition for the DSEs integration and missed its connection with 
the value of the coupling at the renormalization momentum (i.e. the particular value of Aqqd one 
applies to build the solutions) or the critical coupling the scaling behaviour requires to emerge. This 
connection is an important ingredient because it provides us with a manner, through a comparison 
with the physical strong coupling, to discuss whether the scaling critical DSE solution could be allowed 
by the data. 

We will now present, in the following, the comparison performed in ref. [l77l | of the decoupling 
analytical low- momentum expressions, given here by cqs. (|54I59D . and the PT-BFM solutions shown 
to provide a quantitative description of lattice data pot Il87| The main feature in the PT-BFM 
scheme is that it guarantees the transversality of the gluon self-energy order-b y-ord er in the dressed- 
loop expansion, thus leading to a gauge-invariant truncation of the gluon DSE [186]. In this PT-BFM 
scheme for the coupled DSE system, the ghost propagator DSE is the same as the one given by eqs. (l24l . 
where the bare ghost-gluon vertex is approximated by Hi — 1. The gluon DSE is given by 



(1 + R(q 2 )) 2 



9nv 



where 



"l 





(108) 



"2 = 



04 = 



(109) 



In the diagrams of (|109[) for the gluon DSE, Eq. (|108[) . the external giuons are treated, from the point 
of view of Feynman rules, as background fields (these diagrams should be also properly regularized, 
as explained in jl78j |). This is what justifies the four field coupling of two bac kgro und giuons and two 
ghosts leading to the contribut ion ,4. The function 1 + R(q 2 ) is defined in ref. [188| through a so-called 
background quantum ide ntity |178| and can be, in virtue of the ghost propagator DSE, connected to 
the ghost propagator 0; Il87j . The coupled system is to be solved, by numerical integration, with the 
two following boundary conditions as the only required inputs: the zero-momentum value of the gluon 
propagator and that of the coupling at a given perturbative momentum, \i = 10 GcV in this particular 
case, that will be used as the renormalization point. This is done by varying the boundary conditions 
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(In particular, as explained in ref. [I77j . by keeping the zero-momentum value of the gluon propagator 
fixed while cv(/x 2 = 100 GeV 2 ) is ranging from 0.15 to 0.1817) and leaves us with a family of massive or 
decoupling solutions that eqs. (|54I59|) must account for, in the low-momentum domain. The successful 
confrontation can be seen in Fig. 1 of [l77j . from which we extract, as an example, the plot in Fig. 1101 




Fig. 10 Gluon propagators (left) and ghost dressing functions (rig ht) after the numerical in tegra tion of the 
coupled DSE system for a(fi = lOGeV) = 0.16, confronted to Eq. and Eq. If55|l. taken from [177| . The black 
dotted line corresponds to the asymptotical expression including the leading correction and the red dotted to 
the one including the next-to-leading. 



The ghost dressing function at vanishing momentum, F(0,/i 2 ), is also shown to diverge as a — 
a{^ 2 ) — > a CI it, by obeying a power behaviour, s 

F(0) ~ (a crit - a{^ 2 ))-^ , (110) 

where the coefficient k(^ 2 ) is a positive critical exponent (depending presumably on the renormalization 
point, /x 2 ) which governs the tran sition from decoupling (a < a cr it) to the scaling (a = a cr it) solutions. 
This is seen in Fig. 2 of ref. jl77| and confirms the results from the GPDSE analysis in sec. 14.1.41 

Had one let a cr ;t be a free parameter to be fitted by requiring the best linear correlation for 
log[F(0)] in terms of log[a crit — a], one would have obtained a best correlation coefficient of 0.9997 for 
k([i 2 ) = 0.0854(6) and a cr ;t = 0.1822 (which is pretty close to the critical value of the coupling above 
which the coupled DSE system does not converge any more). This last critical value at fi = 10 GeV 
for the coupling can be pretty well translated to that of vIqcd in MS (see for instance eqs. (22,23) of 
ref. [73|]) and the n com pared to in pure Yang-Mills from the lattice B The latter is estimated to 
be 238(19) MeV [Hj, clearly below the former, 434 MeV, for the critical limit for the PT-BFM DSE 
in pure Yang-Mills. 

In summary, one can clearly conclude that the analysis of the coupled DSEs also agrees with the 
existence of both decoupling and scaling classes of solutions and with the pattern for them described 
in sec. 14.1.41 In particular, the analysis of the solutions in the PT-BFM scheme proved the scaling 
one to appear as an end-point for the decoupling family in Landau gauge, when the coupling at the 
renormalization point approaches a critical value. This critical value, at /i = 10 GeV, is well above the 
lattice estimates for the coupling, in total consistence with the conclusion of lattice QCD favouring 
the decoupling solution presented in sec. [3] 

Of course, gluon and ghost propagators and the vertices involving them, altogether with quark 
propagators and the quark-gluon vertex, are basic building blocks to study the QCD bound states. It 

23 It should be noted that the procedures for the lattice determination of A^g mainly work in the UV domain, 
where IR sources of uncertainties as the Gribov ambiguity or volume effects are indeed negligibl e. In fact, there 
are unquenched lattice determinations with Nf — 5 staggered fermions for the strong coupling [1891 ] which are 
pretty consistent with the PDG value. 
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might be that the exact very low-momentum behaviour of gluon and ghost propagators, that we paid 
attention to, is not very relevant for much of the hadron physics. In particular, the quark functions can 
be studied by modelling the product of the dressed gluon propagator and the quark-gluon vertex and, 
regardless of whether t he g l uon is suppressed or massive, the key region for physics is the momentum 
region of p ~ ^Imom f 19li Il92j . Nevertheless, the dressed quark propagator, using massive gluons 
with non-singular interaction in the quark-gluon vertex, do become like expected by the heavy quark 
effective theory (65j . Furthermore, many phenomeno logic al works also appear to support a massive 
gluon solution, as can be seen in the review of r ef. [178| . and references therein, or, as very recent 
examples, in the works of refs. jl93t I194J : Il95t Il96| . This also favours a decoupling solution. 



5 Conclusions 

With this paper, we aimed to give an overview for the current state-of-the-art concerning the in- 
frared properties of pure Yang-Mills QCD Green functions. Very much work has been reported in 
the last few years, modifying essentially the paradigm about this subject and demanding some sort 
of update for past reviews that can be found in literature. About ten years ago, the results from 
Landau gauge DSEs and FRGs analysis agreed with a solution, now dubbed "scaling" , where the 
low-momentum behaviour for the two-point correlators appears to be an enhanced ghost propagator 
and a vanishing gluon propagator at zero- momentum (ctF < and aa > 1, according to the nota- 
tion of Eq. (O). At that time, the lattice estimates for both correlators resulted to be compatible 
with this low-momentum behaviour, other approaches like the ones applying stochastic quantization 
methods and Gribov-Zwanzigcr lagrangian also pointed to the same results and all together matched 
the Kugo-Ojima confinement criterion providing a framework for the infrared solutions of Yang-Mills 
Green functions that was generally accepted. However, more recent results for the two-point corre- 
lators with larger lattices, some of them paying special attention to the problem of Gribov's copies, 
seemed to establish a different pattern for the low-momentum solutions: a finite ghost dressing and a 
finite non-zero gluon propagator at zero- momentum (qf = and ag = 1), which did not agree with 
the scaling behaviour that would require the ghost and the gluon infrared exponents to be related by 
2c*f + ctQ = 0. On the other hand, other authors, applying a particular truncation scheme for the gluon 
propagator DSE which involves some angular approximation for the momentum integration, proposed 
a massive gluon propagator (ag = 1), consistently with the PT gluon propagator. Then, the GPDSE 
was recently re-analysed by exploiting the interplay of DSEs and lattice results and, apart from the 
scaling solution, new ones with up — 0, now dubbed decoupling, were proven to exist, which do not 
observe the scaling behaviour but are totally compatible with lattice results. Thus, both scaling and 
decoupling solutions have been now proven to emerge as solutions of the coupled system of gluon and 
ghost propagators DSEs for different truncation schemes, such as for instance the PT-BFM. The same 
occurs for FRGs. 

On the other hand, we have reviewed the plethora of lattice works on the subject, also discussing 
with some detail the role and impact of the lattice artefacts, and concluded that the current paradigm 
is a massive gluon and a free ghost, i.e. a decoupling low-momentum behaviour for the two-point Green 
functions. Apart from the lattice results, the application of some refinement of the Gribov-Zwanzigcr 
approach and other new approaches also appear now to agree with a decoupling behaviour for the 
low-momentum Green functions solutions. When solving the DSEs, the solutions have been proven to 
be "dialed" by the size of the coupling at the renormalization point, which can be univocally related to 
the zero-momentum value of the renormalized ghost dressing function. A family of decoupling solutions 
corresponds to a family of sub-critical ones for finite values of the zero-momentum ghost dressing and 
for any coupling below a critical value at the renormalization point. As for the scaling solution, it 
can be considered as "critical" , since it emerges for a unique (critical) value of the coupling for which 
the ghost dressing diverges at zero-momentum. A very similar pattern is shown to happen for the 
equal-time spatial gluon propagator and the ghost dressing function in Coulomb gauge. 

Of course, the critical value of the coupling depends on the renormalization point but, once it is 
known for one particular momentum, this value can be propagated to any other by applying the defini- 
tion of the Taylor coupling in Eq. and the scaling solutions for ghost and gluon dressing functions. 
This is of course a consequence of the renormalization scaling for the coupling definition, which does 
not depend on either the cut-off, when written in terms of bare quantities, or the renormalization 
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point, when expressed with renormalized ones. The critical value obtained at 10 GeV in the PT-BFM 
scheme, in agreement with the analysis of the GPDSE with a lattice gluon propagator as input, is 
shown to be definitely above the lattice estimate for the Yang-Mills Taylor coupling, favouring again 
a decoupling solution in Landau gauge. 

The truncation of the tower of DSEs which is necessary to obtain a (finite) tractable system of 
equations implies approximating the vertices. In particular, the ghost-gluon vertex plays a crucial role 
in the analysis which led to find out the decoupling solutions, but also the three-gluon vertex is an 
essential ingredient for the gluon propagator DSE in any scheme. The ghost-gluon vertex benefits of 
Taylor's theorem that has been revised in the appendix, where we also discussed which OPE non- 
perturbative corrections to the vertex the dimension-two gluon condensate induces. For the sake of 
consistency, we also overviewed the results from many lattice investigations about the impact of these 
OPE corrections on the Yang-Mills Green functions and about the possibilities of determining the size 
of the gluon condensate. Concerning the vertices, although some works have been devoted to investigate 
their properties, more lattice and continuum investigations would be very welcome. 

To summarize, in the current state-of-the-art, although both types of solutions are compatible with 
DSEs (and FRGs too), lattice QCD and some continuum approaches, like mainly RGZ, seem to favour 
a decoupling-type of solutions which implies a free ghost and a massive gluon. 
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A A main ingredient: the non-renormalization Taylor's theorem 

A.l What does the Taylor's theorem indeed say? 

A widely used statement, known as the "non-renormalization theorem", claims that, in the Landau gauge, 
the renormalization constant Z\ of the ghost gluon vertex is exactly one. Note that there is no reference to 
a particular renormalisation scheme. Formulated in this way, this claim is wrong. Let us first state and then 
explain below what is true in our opinion : 

1) There is a true and very clear statement which can be extracted from Taylor's paper (the argument is 
given below), ref. [75j |. 

M (-P,0;p) = -z/ p M (111) 

i.e. there is no radiative correction in this particular momentum configuration (with zero momentum of the 
ingoing ghost) 

2) This entails that r^ bc ' Bare (p, k; q) is finite whatever the external momenta, and that therefore Z\ = 1. 

In additi on, w e get also trivially _£ 1 MOMh = x, where MOM;, refers to the configuration of momenta in 
equation (JTTTJ) . In general, in other schemes, there is a finite renormalisation, and this is why one must be very 
careful when using the misleading expression: "non-renormalization". 

3) In particular, one finds in the very extensive calculations of radiative corrections at least two cases of 
MOM schemes where there is a finite renormalisation (and certainly many more) : MOM 9 in the notations of 
ref. [J, and the symmetric MOM scheme. For the latter, we give the proof below. 

The essence of Taylor's argument is actually very simple. In a kinematical situation where the incoming 
ghost momentum is zero, consider any perturbative contribution to the ghost-gluon vertex. Following the ghost 
line in the direction of the flow, the first vertex will be proportional to the outgoing ghost momentum p M , i.e. 
to the gluon momentum — p M . In the Landau gauge this contribution will thus give upon contraction with the 
gluon propagator. Therefore the only contribution to remain is the tree-level one. In other words the bare ghost- 
gluon vertex is shown to be equal to its tree-level value in these kinematics : r^ bc ' Bare (—p,0;p) = — if abc p,i- 
This result has been checked by means of a direct evaluation to three loops in perturbation theory by Chetyrkin. 
In our notations : 

ff 1 (p,0) + ff 2 (p,0) = l (112) 

Note that in the Schwinger-Dyson equation 1)1911 . only Hi is present, and the theorem by Taylor does not tell 
that Hi(p,0) — 1, as seems assumed in many Schwinger-Dyson calculations. 

As an illustration of our point 3), let us quote the formulas from the appendix of ref. Q, reduced to 
the situation we are interested in (£l = 0, rif = 0). The two dressing functions F h (resp. F g ) are defined by 
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Fig. 11 The kinematical situations considered below. The left diagram (0-momentum incoming ghost) corre- 
sponds to /}, below which is known to be equal to one. The right one (0-momentum gluon) corresponds to r g 
and leads to a non-trivial p 2 -dependence 



r^ bc {-p,Q,p) = ~i.f abc r h {p) (resp. r° bc (-p,p,0) = ~if abc r g {p) ) and correspond to the kinematical situations 
depicted in the left (resp. right) part of fig. We have already mentioned that 7h. is exactly one, but this 

does not hold for F g and, indeed, one has at three loops : 
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It is then easy to find the p 2 -dependence : 
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In ref. [6(| the non-renormalization theorem is understood as the statement that the vertex reduces to its 
tree-level form at all symmetric-momenta points in a symmetric subtraction scheme. How ever this statement 
is not supported by a direct evaluation. Using the one-loop results of Davydychev (ref. 78]) one gets in a 
symmetric configuration the value 



r^{p,k; q )\ P ^= q ^=-it bc {p» (l+ £^0+^))+^ ^^0+^)} (H4) 

with 0= ^C7 2 (f),C7 2 (f) = 1.049 

According to ref. [6fJ] the coefficient of p M should be one. The presence of a 3 in the above formulas implies on 
the contrary that the vertex will in general depend on the momenta : using the results given in the appendices 
of ref. [J one finds for the leading p 2 -dependence 

-'/"ifi'£> ,1 °4 , ( (9+ l* , '" +(3+ i*»)}' 

This dependence is logarithmic, as is expected in a perturbative approach. Furthermore, in ref. ( 66]) it is 
supposed that the vertex function takes the form (q 2 ) e (k 2 ) m ((q — k) 2 ) n with the restriction £ + m + n = 0. This 
restriction comes from the assumption that the symmetric vertex is equal to 1 for any p 2 , which, as we have 
just seen, is actually not the case. Therefore, it should be necessary to adopt a more general point of view and 
keep open the possibility of a non perturbative effect on Hi. We should mention that, actually, the problem 
of the p 2 -dependence of the ghost-gluon vertex has already been addressed in refs. [1821 : 1 183 ] . However these 
aut hors work under the condition 2of + Qg = which appears not to be satisfied by lattice dat a. Also in 
ref. [1651] . the impact of the OPE non-perturbative corrections, as the one we dealt with in sec. 13.81 is studied 
for the ghost-gluon vertex in order to go beyond the approximation of taking the Taylor kinematics for the 
transverse form factor. We will briefly pay attention to this in the next subsection. 



A. 2 Non-perturbative corrections for the ghost-gluon vertex 

The non-perturbative effect resulting from a non-zero {A 2 } in a OPE expansion, shown in sec. 13.81 to have 
non-negligeable effects on ghost and gluon propagators at energies of the order of 2-3 GeV and still visible 
at a roun d 7 GeV, can be also advocated to have an impact on the ghost-gluon form factors introduce d in 
Eq. ()18p . in particular on the transverse form factor, Hi, needed for the integration of the GPDSE in sec. 14.11 



46 



The procedure outlined in sec . l3781 is also in order to compute the Wilson coefficients for such a non-perturbative 
contributions. Then, Eq. (|86|l can be particularized to be 



Hi{q,k) = Hr\q,k) (l+ ^tT^ <^V)>ms) 



(115) 



q-k 



-V 




(a) 



(b) 



Fig. 12 Diagrams contributing (altoget her w ith their appropriate permutations) to the ghost-gluon form 
factors for the proper ghost-gluon vertex |165lj . 



where, although it is not divergent, one can apply a finite renormalization to require Hi to take its tree- level 
value, Hi = 1, at a given momentum scale, fj 2 , for a particular kinematical configuration of p and q (for 
instance, p 2 = q 2 — fj 2 ). This dependence on fj, s houl d be understood for the Wilson coefficient and t he f orm 
factors. The Wilson coefficient can be obtained [1651] at tree-level by evaluating the diagrams in Fig. [12] and 
reads: 



Hf el \q,k) 64 M ' V q 2 (q-k) 2 fc 2 (g-fc) 2 

Thus, after modeling the ghost-gluon form factor all over the range of their moment a by t he insertion of some 
infrared mass scale to saturate the powers of momenta in the denominators of Eqs. (|116[) (as a simple way to 
avoid the non-physical divergence coming from the se inve rse powers), one obtains the model for the ghost-gluon 
transverse form factor that was mentioned in sec. 14.1.21 This model provides us with a first correction to the 
usual hypothesis for the GPDSE integration: Hi = 1 that is encoded by the OPE expansion through (A 2 ). 

It is interesting to notice that, had we considered the kinematic configuration for the Taylor scheme, k = 0, 
one would obtain a non-vanishing OP E po wer correction for the transverse form factor; although no deviation 
from the Taylor result shown by Eq. (|112[) would result, because 

Hi(q,0) + H 2 (q,0) = Hr t (q,0) + Hr\q,0) = 1, (117) 

as it is proved in ref. [165lj . 



B The Dyson-Schwinger equation as a Ward-Slavnov- Taylor identity 

A very general method to derive Slavnov- Taylor identities consists in taking advantage of the transformation 
properties of 



V(A) det A^exp 



d x £ 



2a 



(118) 



under gauge transformations (cf. |76lj ) . 

M. is the Faddeev-Popov operator and the notation <, >j indicates that the source term J has to be kept, 
although it will eventually be set to (this is denoted in th e following by the suppression of the J subscript). 
Taking the derivative of the gauge transformed of eq. (|118l) with respect to the gauge parameters leads to the 
general Slavnov- Taylor equation 



(y)F ( - 2)ba (y,x) >j 



(119) 
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F^ ba (y, x) is the ghost propagator and its presence here is simply due to its very definition as the inverse of 
the Faddeev- Popov operator. If one derives eq. (|119p with respect to J p (z) one gets : 

1 

a 



- < (d p A d (x))A d (z) >j = < D db (z)F {2>ba (z,x) >j 

+ < f d 4 yj^(y) D cb (y)F^ bd (y,x)A d p (z) >.j 



(120) 

A first consequence of this relation is the triviality of the longitudinal gluon propagator. To see this, it suffices 
to derive both its sides with respect to z p and to set J to zero. The result is 

i < (d^A d (x))(8 p A d p (z)) > = < d p D db (z)F < - 2)ba (z y x) > 

= 8 ad 5 i (z-x). (121) 

In order to derive the second line we have invoked the fact that d p D p (z), the Faddeev-Popov operator, is the 

inverse of the ghost propagator F 1 ' 2 '. Thus, in momentum space, the general form of the gluon propagator for 
an arbitrary covariant gauge reads 



G$ a "(q) = 5 a 



G {2 \q 2 )[8 pv - 



q P q v 



Turning now back to eq. (|120[) and letting J go to zero we obtain 

- < (d p A;(x))A d p (z) >=< Df{z)F^ ba {z,x) > 



(122) 



(123) 



which is nothing else than the GPDSE. Actually its l.h.s. involves only the longitudinal part of the gluon 
propagator, that we have just seen to be trivial : 



- < (8 p A;(x))A d p (z) > d p D-\x,z) 



(124) 



where the □ symbol stands as usual for the d'Alembertian. As for the r.h.s it can be rewritten as : 

< D db (z)F {2)ba (z,x) >=< d p F (2)da {z,x) > +i < gf deb A e p (z)F {2)ba (z,x) > . (125) 

The 3-point gluon-ghost Green's function can be expressed in terms of vertex functions and propagators 
through 



G 



('i)fgh 



(p,q,r) =-if d 4 x d 4 t d 4 ze lpx e lrz e ,qt < A f p (t)F (2)Bh (z,x) > 



(126) 
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f f9h r„(p,r;q)(2Tv) 4 84p + q + r) 



Now, we Fourier-transform Eq. (|123|) , use eqs. (|124|126p and obtain 



^F(k 2 )~gf dsb f 
G(q 2 



bf eba f d 4 q F(k 2 )F((k + q) 2 ) 



(27T) 
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(127) 



where the usual form of GPDSE can be recovered from by multiplying with k p and dividing by F(k 2 ), which 
leads to 
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r v (k,-k ~ q;q) 



(128) 



This is a general result, valid in any covariant gauge. Of course the a-depending (longitudinal) term disappears 
in Landau gauge. r v (k, —k — q; q) is related to the quantity, F^, that was previously introduced in section [T] 
through 

r v (k, -k — q;q) = —igky.r^k, -k - q; q) 

and it is usually decomposed into F v (k, —k — q;q) — g [k v Hi(k, q) + q v H2(k, q)]. After inserting this in eq. (|128jl 
and restricting to the Landau gauge case one finally obtains 



F _1 (fc 2 ) = l + £/ 2 iV, 



r d 4 , 

J (27T 
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(129) 
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